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ABSTRACT 


This thesis is the result of literary research on the subject of 


sensitivity in linear feedback control systems. It is a synopsis of 


information obtained from various technical publications and is designed 


to give the reader information on the theory of sensitivity and its 


application in design problems. Included are 


(1) Definition of sensitivity function and root sensitivity 


(2) 


(3) 


for incremental variations of single parameter and their 
design application. 

Definition of sensitivity function applied to the large 
parameter variations and its design application. 

Special topics such as specification for sensitivity 
function, sensitivity integrals, and introduction to 


; # sy 
multi-parameter sensitivity. 
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I 
INTRODUCTION 

One of the foremost properties of feedback is its ability to 
reduce the sensitivity of a system to variations of the system para- 
meters. Suppose there is an element in the system which must be used 
because of its special properties. However the element is sensitive 
to environmental conditions (such as temperatire and pressure), or is 
subjected to aging, or has wide manufacturing tolerances such that when 
it is replaced the pszrameters of the new element may be markedly diff- 
erent from those of the old element. Whatever the reason may be, 
suppose the values of the element parameters vary significantly. It 
is desired that despite these variations certain system properties, 
such as the input-output transfer function or output impedance should 
remain substantially constant. The troublesome element could perhaps 
be replaced by a better one, but this may be very expensive or imposs- 
ible. In any case we would like at least to liave an alternative to 
replacing the element. In most cases the troublesome elements are the 
active elements, vacuum tubes, transistors, energy conversion devices 
such as motors and generators. ut they may also be passive elements 
or transducers. In general the system assumed in an analytical study 
will never precisely match the actual physical system; it is therefore 
important to know the effects of possible variations. 

A general knowledse of the effects of parameter changes can te 
used as a guide to system modifications which would improve the over- 
all performance. For all these reasons measures of closed loop system 
sensitivity to open loop system parameter variations are an integral 


part of the analysis protiem. 


Notions about system sensitivity were in the forefront when the 
feedback concept was initially developed. This was natural, even 
unavoidable, since feedback systems possess the fundamental physical 
property that the effects of variation in the forward loop, whether they 
are taken as changes in open loop transfer function G(s), or as depart- 
ures from strict linearity or from freedom from extraneous noise, are 
reduced by the factor 1 + G in comparison with the effects which would 
be observed in a non feedback system. Accordingly, sensitivity measures 
were indispensable to any rational discussion of feedback systems and 
a useful, classical definition of sensitivity was made. Except for a 
minor modification the definition of sensitivity, given by Bode, remained 
unchanged for over a decade. Perhaps this static nature surrounded by 
dynamic growth in most other areas of feedback systems engineering, 
made the concept fade with little exposure to sensitivity concepts. 

In recent years classic sensitivity has become a more popular subject 
in automatic control and also in network synthesis. Finally, the 
emphasis on pole-zero specifications for system characteristics gave 
rise to new concepts of sensitivity, with associated new measures. 
These measures called here "gain," "open-loop pole," and “open loop 
zero" sensitivities, were evolved to account for changes in the posi- 
tion of closed loop poles due to shifts or changes in open loop gain, 
poles or zero. 

This thesis is the result of a literary research on the subject 
of sensitivity in linear feedback control systems. In section II, the 
system sensitivity and root sensitivity, for a system subjected to 
incremental parameter variation, are discussed. Based on the meaning 


of open loop pole and zero sensitivity, a design technique for 


10 


compensating the system is also introduced. In Section III, the large 
parameter variation and a new definition of sensitivity is presented. 
It is shown that when the leakage transmission between input end output 
is negligible, the new definition and the classical definition of the 
system sensitivity are identical. In Section III, an attempt is made 
to utilize Nichols Charts for evaluation of the sensitivity at a part- 
icular frequency range. A design philosophy and technique is also dis- 
cussed. Section IV discusses some special topics which are of parti- 
cular interest today and could be selected as interesting topics for 


further investigation. 
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II 
SENSITIVITY AND INCREMENTAL PARAMETER CHANGES 
II-1 - System Sensitivity and Root Sensitivity 
Several definitions of sensitivity have been offered in the past. 
Bode defined the sensitivity of the overall transfer function T to the 


plant parameter k as 








2k 
st. tk , k 
k 3 Ln T AT 
T 
Horowitz used the inverse of Bode!s definition, 
oT 
gt , AT y T 
k § Ly k ok 
k 


Defined one way or another, sy is generally known as "classical sensi- 
tivity" or "system sensitivity," since it involves the system transfer 
function T. Several other definitions have also been suggested, and 
each one of these definitions concentrates on some particular engineer- 
ing or mathematical simplicity. A. Schulke suggests since sensitivity 
relates changes in the transfer function of a system with respect to 


changes in parameter of the system, then one can define, 


ee | eee 
sensitivity S >) bn x 


and in general T is a complex and may be expressed as 

T(x) = A(x) + jB(x) 
where A(x) is the system attenuation function and B(x) is the system 
phase function, and from that, 


$= HR + ae {ae + sven} 


x 





= Xx AGO op jx —O3QO_ 


ax 
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Then assuming that x is real, it is apparent that the real part of 

Z is directly related to the system attenuation and that the imaginary 
part is directly related to the system phase function. The use of 
s@ which Schuike calls it “sensitiveness" permits comparison of both 
attenuation and phase functions relative to changes in a common parea- 
meter. 


The most adcopted definition of system sensitivity these days is 


the one sugsested by Horowitz, i. e. 


gt ag ene 
k din k 


and we adopt it through the coming discussions. 

In recent years with the increased utilization of the pole-zero 
approach, it has become increasingly important to examine variations 
in position of the poles and zeros of the network function due to 
changes in network parameters. This is known as the "root sensitivity." 
Formal definitions of root sensitivity vary from author to esuthor. 


Some define it as, 


= ’ 
Se eee! 
x ax ox 
x 


where (-s,) is the root of the characteristic equation and x a para- 


j 


meter which could be the gain K, or a pole or a zero. Huang uses the 





definition 
3 
8 — = ds, 
Sy a 3 ax 
Ox j 
x 


Mc Ruer ard Stapleford prefer different definitions for sensitivity to 
s § ¢ 


gain, dh and sersitivity t. poles, Sa or sensitivity to zeros, So : 


Again each one of these definitions leads to a particular engineering 
or mathematical simplicity utilized by that particular author. 
II-2 - Sensitivity Relation in Control Feedback Systems 

It has been noted by Bode that for any linear system the transfer 
function T(s,K), relating a response to an input function can be in 


the bilinear form of, 





C(s) (s) 
aT 
and sr . 7m : es 
2K 
K 


Both numerator and denominator of T are of the form, 
H(s,K) = q(s) + Kp(s) 

where K is a variable parameter and p(s) and q(s) are polynomials. The 
roots of H(s,K) = 0, which characterize H, depend on the value of K. 
H(s,K) = O constitutes an implicit relationship between s and K, and 
as K takes on all positive values, a number of curves in the S-plane 
are obtained. These curves are the loci of the roots of H called the 
"root loci." . 

The root locus was introduced by Evans in 1948. He used it main- 
ly for the investigation of control systems with a transfer function 
T, in the form of, 


G = Kp(s 
1+c¢ q(s) + Kp(s) 


T = 
where G = Kp/q is called the forward gain in a unity feedback con- 
trol system. For such systems, only the denominator is of interest. 
As the result, there arose the generally accepted terminology in which 


the zeros of p and q were called open loop "poles" and "zeros," 


respectively and those of 1+ G = 0 are called "closed loop poles." 
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II-3 - Relationship Between System Sensitivity and Pole Zero 
Sensitivities 


The transfer function T(s) can generally be written in the form 


of 
n 
SiM(s +2 ) 


T(s) Mw an (2-1) 


TT(s + 8,) 
1 j 


and let Z s 
i) Aut 


ature, gain of some element, etc.). Then 


and g be functions of some parameter x (such as temper- 














n m 
fnT = ng + Bh(s + z,) - Sls + 8) 
1 1 
and differentiating above we get 
n m 
ats de} hag+ 3 24 Saad 
T gx 1x st 2 1 dx s+ . 
dividing both sides by dx 
x 
= = Dh. Se Se 
a Ye Ss ->ax 4 1 (2-2) 
= a= l ox 8 + 2, 1 dx st 8 


Equation (2-2) gives the relationship between the classical definition 
of sensitivity and the closed loop pole and zero sensitivity, where 


the closed loop pole and zero sersitivities are defined as 


& 
‘ 2 
So = x . (2-3) 
ax 
ff Fane 2) (2-4) 
Ls ae 
ox 


This definition is not necessarily the most useful definition at ail 
situations. 
Pole Sensitivity (closed lioop pole): 


The root locus is defined by: 
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D 


H(#,K) = q(s) + Kp(s) = O (2-5) 
hence: 
n 
IT (8 + py) 
R= a) = - tL _ (2-6) 
p(s) TI (8 + 2.) 
i=] 
and 
n m 
Bn (K} = 2 bn(s + P,) : 2: Ln (8 + 2Z,) + 15 (2-7) 
j=l i=l 
on dK n 1 m 1 
*K 7 $12 Gap) - 2 Grey (2-8) 


K = - q(s) 
p(s) 
then 
dK . __a'p- p'd (2-9) 
ds 2 
p 
and 
UP = p'q_ 
dK 2 ; : (2-10) 
Co a ae ine 
x q q p 


For simplicity of notation Equations (2-6) and (2-8) can be written as, 


K = -TW(et+s er (2-11) 
i 
che yet ds yy a (2oheo 
K iis+s 
i 
which has to be evaluated at s = “Bae The roots of H = O are 


denoted as (8) to distinguish them from the variable s which may of 


course assume any value. In other words -s, is ao point on the root locus. 


j 
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And 


een 1 if ~ S. is a root of aq(s). 


i = cl if - Ss, is a root of p(s). 


From Equations (2-10) and (2-12) we can write, 





3 
Sia = nat = l = “pq 
dK oS aay ap’ = pq' 
‘ s+ s 
‘ S=8 
j 
zs > K 
q' + Kp‘ (2-13) 
s 
where s.7 is the closed loop pole sensitivity with respect to K and 


sometimes referred to as just pole sensitivity with respect to K. 


Furthermore the residue of a function A(s) at pole s*a is defired as 
B(s) 


Residue _A(s)__ = lim As = a) A(s) 


at s=a B(s) sa B(s) 


which is true for simple poles. Noting that 3H = q'‘' + Kp‘ and 


is 
applying above residue theorem we get : 
z 
Se = - _ Kp os © Xp = Q, (2-14) 
q' + Kp! 2H 
o8 


where Q is the residue of 


Kp(s 
T(s) = G(s) + Kp(s) 


at point s = *,. The same result can be obtained directly from (2-5) 


<3 ds aE 

Se = j = oo» K 3K (2-15) 
dK Te 
K a84 


and since H = q + Kp, then 


3H 
- on 
oJ = -K_K = - K = - Kp = -2, 
K qH dq + K dp YH 
Xs, ds ds Is 
sm-8) e=8), 


The relation 3 & -Kpiq' + Kp’) is quite convenient since it oniy 
requires knowledge of the "open loop" transfer function and need not 
be in factored form. The evaluation of the pole sensitivity 53 by 
the residue of the closed loop transfer function will call for the 
complete determination of the closed loop transfer function. 
II-4 - Slope of Root Locus 

Slope of the root locus is the tangent of the argument of 5 


and since, 


rs b 
s,) Ee 8 (2-16) 
AK 
K 
°5 
Then ds. = SS.” AK (2-17) 
j K K 


Since AK is a positive number, then 


K 
VA Le (2-18) 





where /As, represents arg ds,. So 
dy i (oe oe 
qx ™@ tan/ds, = ten Ss, a s (2-19) 
Res,,j 


(2-20) 
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Input Output 


Sp G(s) 


Fig. 2-1: Unity Feedback System. 





Fig. 2-2: Root Locus Slope. 
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If we write 


- je3 
s+ Pp, 4, e 








i 
and since 
Ii si - I, i: ws sin © 
5 5 Z 5730 a 
+ 
oP j j 
then m sin 6 n sin 8, 
Le Be 
dy = =] i=] i (2-21) 
dx m cos @ n cos @ 


Equation (2-21) could be written in simple fom, 


sin 9, 
dA, i 


dy = 


dx vA cos Q, ree) 
2 q 
where A, = 1 if 4, and 8, are due to poles 
Ay = -l ig {, and e, are due to zeros, 
and Me = magnitude of ch ¢, 
eo. = argument of soa 8; 
Then 
» iti Vy 
DAA & - x)? + (y - y,)? 
dy = i i i (2-23) 
dx x-x 
fe) ee 
2 2 
Zul -x +0 - 7) 


In some cases it is easier to compute dy/dx from the slope of ds/dK 
obtained by differentiation of - K = p/q as shown in example discussed 


at the end of this chapter. 


20 


j 
*) 
l 
a.~ 
— 
| 
| 
| 
! 
s), 8, ! 
] 
| 
| 
e 
8, 2 
Fig. 2-3: Argument end magnitude of 7 we 





Fig. 2-4: (a) Definition of K 


8 
(b) Evaluation of 
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8 
II-5 - Graphical Evaluation of ee 


The pole sensitivity at any point "A" corresponding to a particular 
value of s can be evaluated graphically by summation of angles and withe- 
out evaluation of K. Choose a point C such that AC is perpendicular 
to the locus and evaluate the phasor AC = AS. Evaluate the value of 
the phase of K at C by summing angles e, and let it be 7+4a. This 
procedure is based upon the independence of the derivative of an 
analytic function on the direction of differentiation. Hence AK is 
chosen to be A'C' which is mapped into AC. 


Then: K # j(Act)K 





gs = ss = As = 8, (2-24) 
jos §QK 0 AK 
jk 
“3 
Therefore knowing the pahsor As and the angle 4%one can evaluate sy : 
II-6 - Curvature, Pole Sensitvity and Their Relation 
The curvature at any point on the root locus is proved to be 


(Appendix II), 


1 2. - im'"*),» = = t0\“").» (2-25) 
Sam Imj(tK)_ Re(4aK) | 


At first glance it seems that the curvature and sensitivity are someliow 
related. An example, discussed at the end of this chapter, will reveal, 
however that the curvature in that example is constant while the pole 
sensitivity is equal to cot $ and consequently is a continuous variable. 
The independent nature of the curvature and pole sensitivity can be 
understood by noting Equation (2-25) and the definition of sensitivity. 
Pole sensitivity depends on first derivatives, while the curvature 
depends on the second derivatives, and in general the two derivatives 


are independent. As suggested by Ur we can think of curvature as 
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"kinematic" property and the pole sensitivity as a “dynamic 
property." 
II-7 - Open Loop Pole and Zero Sensitivity (Figure 2-5) 

Let G(s) be the forward transfer function and L be the loop 


transfer function. Then L = GF and the system characteristic equa- 


tion is 
l1+L = QO 
and if 7 is the system root, then 
1 + L(s) = 0 
ons’ 
where 8 = -s, represents the point on the root locus. If K is the 


j 


gain constant of L(s) and Z,, P, its open loop seros and poles, then 


1’ =i 
we can write: 
L L(s,K,Z,,P,) 


and taking the total differential of L: 


ni m 
dL = Wie + Wy + ZF wb. + Se : 
ys ox i=l 2Z, dz, i=1 2P, dP (2-26) 


On the root locus 1+ L# 0 or L(s) ® -1 and therefore L(s) = constant, 


and the total differential dL is zero for 5 @ -s Then letting 


j" 
dL = O and 6 = fi : Equation (2-26) gives, 
0 = wk dam > dz, +) dP 
38 | © Se MZ i 
ans, Seog i <a i ems, 














Rearranging we get 


"5 (38) anes Piss, 1) bs +504) 2 a DG ibs, ae 


But bs, itself is a function of K, Zs and Pi so 


a = eae Zi. P) (2-28) 
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Taking the total differential of e we vet 
K Xs 8 ds 
ees ON S'  fedZ es dP : 
ds, eet tm tt fw (2-29) 
Equation (2-29) suggests that oe be written as 
= s 2; 
ds, = S! ak + S's)? az + s J ap (2-30) 
j K = Z i P i 
K i i 
Equations (2-29) and (2-30) thus define the sensitivity of root a 
to gain K, as 
*; D Ys ‘ 
Sy = j (2-31) 
oK 
K 
which is the same as Equation (2-13) and the sensitivity of root S; 
to open loop zero Z; as 
s 
s,) 8 os (2-32) 
i 62, 
and the sensitivity of root 2 to open loop pole Py as 
8 er - 
ges ee (2-33) 
. P 
{ 
S S s, 
II-8 - Relationship between S,°, 8,7, and oo 
Comparing Equations (2-27) and (2-30) yields: 
4 
Ss K (217K) (2-34) 
oL/3S J 8 = a3 
6 
si = (ee) (2-35) 
i L/dS S$ = a 
s 
s,) - ae) (2-36) 
i L/S g = -8, 
n 
ee TT (s + Z,) 
L #*= Kl (2-37) 
m 
T (s + P,) 
1 











Then n 
TT (s + Z,) 
+4 Pec es, (2-38) 
¥K “94 m (+) aS 
TT (s + P,) j 
l ;s = -8 
j 
But since 1 + L(s) = 0, 
8 = <§ 
j 
then L(s) = 1, and Equation (2-38) yields to 
s = ==) 
YL = -1 (2-39) 
YK/s sB ax K 
Then (2-34) becomes 
8 
S — Be//\1. (2-40) 
; Ys) 8 = -s 
j 
Similarly ne-l 
KIT (8 + Zz.) 
bs 
3Z m 
ij s = -s 
j a (s + PD) Se 
= L = -l (2241) 
s + Zs ea “8, as + Z; 
Then, 
es Lista ey 7 1/(-B + A 
L/AS & = se: QL/dS) . ae 
8 
oS sj = 5xd (2-42) 
das + Z; 
Similarly 
s as 
s,) = K (2-43) 
i 24 = Py 


Equations (2-42) and (2-43) dictate that the sensitivities to all 


8. 
singularities are directly proportional to e and inversely pro- 


portional to the distance between -s, and the singularity involved 


j 
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(in this case s = -Z, and S = “P. are the signularities) Thus what- 


8 
7 j 
K P 
A particular case of Equation (2-43) is for Pa = 0, i.e. the 


1 


s 
ever properties are found for S&S J may be extended to s,; or S$ 


root sensitivity to the pole at origin of the S-plan. 


8 
‘, 5 5 
S = K (2-44) 
Fo (s,) 
j 
< 8 
Then So is proportional to Sy , the constant of proportionality being 
Oo 
1 , a complex quantity. 
8 
j 
II-9 - Relationship between Root Sensitivity of es and Residue 
In Equation (2-14), we proved that in case of unity feedback, 
8 
IPSs = 
SK Q, (2-45) 
where Q, ig the residue of T(s) = Kp(s) at point ra 
q(s) + Kp(s) 
In this section it will be shown that if Al is a single root, then 
8 
j 
S = F(-s 2-46 
P (-8,) Q (2-46) 


where F(-s,) is the feedback transfer function F evaluated at s = -s 


j 


and Q; is the residue of the system transfer function at Ce) For 


j 


unity feedback and simple root F(-s 5) = -1 and Equation (2-46) yields 
to Equation (2-45). The overall transfer function T(s) of the system 
shown in Figure (2-5) is 


G(s F(s)G(8s 


Ts) = TY ¥F(s)e(s) ~ “FCs)tl + F(s)C(s)] (am 
But L(s) = F(s)G(s), then 
COR oe) ae (2-48) 


F(s){1 + L(s)} 


The residue of T(s) at -s, is 


j 


Qy = (8+ 8,)7(8) all aha 


a= -s, ~ “R(s)fl + L(s)}, . “8, 
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Then, let the rightmost expression be denoted as Ry» i.e. by 
definition, 
R ,(-8 ss 2-49 
rh > Q, ( ) 
Then RAFCL + L) = (8 + 3 ,)L 


Taking the derivative with respect to s of both sides, 


aE (1+ LR, + FR, + FfL+ Ly °%y - L+ (s+ 9) $5 


8 
At ss “8, i.e. at a point on the root locus, L = -l the above 


equation reduces to 


wL 
Biaor aa m 3] 
8 = -8 
OR R,(-s,) © l/F(-s.) | \L (2-50) 
j  § j | oxi 


Comparing Equations (2-50) and (2-40) we get 
8 
J = F(- 2 
Sy F( 84) Q, (2-51) 


Now turning to the case of N order root at -8, we can expand 


j 


L 
ne 


into partial fraction as 


Q Q Q 
T(s) = aa Ae + —atsjr tt 


$+ 8 


j (8 + 8;) 
+ terms from other roots (2252) 
where 
(N = K)!? \e F(1 + L) 
8 & -§8 
j 
Again defining the quantity inside the small bracket as Ry» then 


repeating the operations as before, one obtains after repeated 


differentiation 
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- N! 
oF SEU EENRGIES Ra SRE (2-54) 





F(-s,) L 
oN 
8 8 on 
and from Appendix IIT 
a (-1)" WN! (2-55) 


\s* | s = “8, 


Theorem 
When at is a single order system root, the sum of the 
sensitivities of EF to open loop poles and zeros is equal to unity, 
i.e. 
Ys! + 53 = 1 (2-56) 
i i i i 


This is easily seen by referring to the construction of root loci. 
If all open loop zeros and poles are displaced by the same amount 
2, then all closed loop roots are displaced by the same amount a, 
that is if dz, = dP, =§ for all i, then ds, «4 for all i. There- 


2 J 


fore if in Equation (2-29) we put aS = dz, im qgP, = and dK = 0 we get 


ds is 
_ + 5 = a 
i i 
A more rigorous proof is as follows. Rewriting equations (2-13) 


and (2-42) and (2-43) we get 


s S, 

i © Sx 
i Z - 84 

8 Sj 

ae = Sx 
i - = P. 
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s | 


The above relation is only valid for simple poles. 
II-10 - Graphical Method for Determining Open Loop Pole and 
Zero Sensitivities 

Equations (2-42) and (2-43) suggest that to singularities which 
are close to system root (-8,), as is more sensitive and for singular- 
ities which are far away, S5 is less sensitive until it becomes 
insensitive to singularities at infinity. This concerns the magni- 
tude of sensitivities. But sensitivities are vector quantities, 
since the change of a parameter may move the roots in different 
directions. It is then helpful to make use of Equation (2-56) toget- 
her with Equations (2°42) and (2-43). 

Referring to Equations (2-42) and (2-43) where Z and P indicate 


open loop singularities and -s, is the system root in question, we 


J 


can draw a vector from “a toward each and every pole, and away from 
each and every zero. The length of each vector is inversely propor- 
tional to the distance from = to the singularity concerned. Then 
construct the sum U of all these vectors, which is a vector itself, 
(Figure 2-5). If U is taken as a unity vector in magnitude and phase, 


then the other vectors measure the sensitivity of of to each 
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Fig. 2-5: (a) Construction of vector diagram 


(b) Vector addition to obtain U. 
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Fig. 2-6: A Basic Feedback System. 





Fig. 2-7: Phase Relationship. . 


31 


singularity respectively. So we regard it as 
u =1fo = 1° (2-57) 


and is named by Thaler and Rung as "unit-sensitivity" vector. This 
8 
unity scale does not apply to sensitivity to gain, ee 


this technique is quite similar to the graphical method of evaluating 


Actually 


the residues for simple singularities, and by putting U = 1 fo we 
establish the necessary reference scale. Phase of sensitivity vector 
is measured as positive in the clockwise sense starting from the U 
vector. This seemingly arbitrary sign convention in fact comes 


from Equations (2-42) and (2-43) which are the basis of this vector 


diagram. A “ 4 
8, s_J s_J s J 
Sp = K = K ws K (2-58) 
i 5, o P, (-P, *) cA 


The denominator is the vector from root (-s,) toward pole (-P,) as 


j 


shown in Figure 2-7. The phase relationship of the above equation is 


8 
j 8 
$ & j @ —_—> 
(es es: [1 (2-59) 


Relation (2-59) is true no matter what conventions are applied to the 
measurement of the angles. Since the U vector has been found to be 


) 3 8 

equal to 1 /o as far as S J and § J are concerned, phase of S J 
Pa 24 Ps 

will be measured starting from U as zero phase. On the other hand 


—> 


and vA are measured in the conventional way, i.e. 


8 
j 

Sx 
starting from the positive real axis and counting positively 


counterclockwise. . 


j 
In Equation (2-59) the quantity Sx does not depend 


on i, i.e. it is the same for all i's. Thus for each i, 


8 
S j 
Ps = constant - V; 
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ptr that the larger the value of / Ei the smaller must be 
= 
Spy . Since is measured conventionally (positive 






counterclockwise) then 
wise. 


S, . 
Next we prove that vector $7 lies on U. From Equation (2-44) 





Ss, 
oe 7 j 

ie 
~ re Li ae 
Figure 2-7 shows the angles ie > S * =QA , the latter 

Po Ss, 
S 

being measured from U. Since as shown in Figure 2-/7, Po “+r 


[s\ gives the direction of U on the S-plane, it is thus establish- 
Ss, Ss. 
ed that on the S-plane, 5S J always lies on U. Since Se indicates 


K 
the direction in which the root moves when K varies, i.e. the direc- 
tion of the root locus, the above result can be stated as: "At any 


point on the root locus, the U vector is tangent to the root locus." 


Also Equation (2-59) reduces to 
s J 
Ps = - V5 (2-60) 


II-11 - Locus of U on the S-Plane 
Consider the plant to be compensated, 


Z 
a K(s +1) . 
ee s(s + Pi )(s + P,) Gao) 


and the desired dominant-root location = as indicated in Figure 
2-8. A spirule measurement shows that an additional phase of + 9 
is needed at location a Thus a lead network with a zero at Z and 


a pole at P is needed. The question is, how does the tip of the U 
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vector move on the S-plane, when Z and P take all possible values 
on the negative real axis? 

Now sensitivity vectors should be drawn for all the plants poles 
and zeros, ignoring for the moment, the need for compensation. All 
these vectors can then be added to obtain the vector sum QI, which 
may be called the "uncompensated" unity vector. Then the geometric 
locus of U is a circle with a center I and a radius of r = 1 : 
This circle is named the (U) circle. 

QI would be the unity vector without "P" and "Z" (Figure 2-8a). 
When P and Z are added so that POZ = b, a sensitivity vector must 
be assigned to each of these singularities, Qm= 1 is assigned to 


QP 


the pole and Qu= 1 is assigned to the zero. The vector QI then 
increases by the ie quantities IM = Qm and MU = Qu. QU is the 
final U vector. Note that, in Figure 2-8a, IN and MU are two equal 
but opposite vectors. 

Figure 2-8b shows that when P moves along the real axis, since 
Qm = 1 , m moves on a radius of r = 1, because of inversion of 

QP ad 
real axis with Q as center and ratio of inversion one. Such a circle 
is referred to as the (m) circle. If we translate this (m) circle 


by QI, then we get (M) circle which is the locus of point M. So the 


locus of point M is the (M) circle, or radius R = _1 and with I as 


2d 
its upper most point. Figure 2-8c shows that: 
“a 
IU = MN = 2R sin MON 
2 
IU = 2R sin f = 1 sin = constant (2-62) 
d 


Hence U moves on a circle centered at I and of radius r= 1 sin >. 


So the following result can be stated: 
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Fig. 2-8: Geometrical Locus of 
Point U. 
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"The geometric locus of the tip of the unity vector is a circle, 
centered at I and of radius r = 1] sin b, where I is the tip of the 
‘uncompensated’ unity vector, d 7 imaginary part of the dominant 
system root and @ the phase shift to be introduced by the compensation." 

Note that there are limits to the geometric locus of U on the 
(U) circle. One limit corresponds to the extreme case where P = om, 


the other is for Z= 0. Later it will be shown how to obtain these 


limit points on the (U) circle. 


II-12 - "Use of the U Locus" 

Starting with the uncompensated plant and a desired location for 
the dominant root, it is always possible to be able to draw the (U) 
circle and be able to obtain the two limit points on the U locus on 
the circle. 

Conversely after a point U has been chosen on the U locus, it 
is possible to perform the graphical construction in reverse order 
and thus drive the corresponding Z and P, i.e., obtain the compensation 
needed, The following will explain how to choose the U point on the 
U locus, which is the seence of this construction technique. 

In Figure 28a, note that OJ is perpendicular to IU where J is 
the midpoint of arc MN, P and Z can be obtained by drawing OJ per- 
pendicular to IU, cutting the (M) circle at J, drawing arcs JM = JN = 
®. P is then determined by drawing OP parallel to IM, and Z is det- 
ermined by drawing QZ parallel to IN. 

An example is given in Figure 2-9a, where a desired U point is 
shown (an arbitrary choice for illustrative purposes). OJ is drawn 


perpendicular to IU. Next, are JM = arc JN = p is drawn. In Figure 
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U circle 


(a) 








Fig. 2-9: Determination of P and Z for given U. 
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(a) 


(b) 


Fig. 2-10: 





U circle 
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Determination of the limit points of geometric 
locus of U. 
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2-9b P is obtained by drawing QP parallel to IM and Z is obtained by 
drawing QZ parallel to IN. 

Now the problem recurs of determining the limit points of the 
U locus, using the same construction procedure. The extreme right 
limit point corresponds to the case where one compensator's singu- 
larity (here, the zero) is at the origin. The extreme left limit 
point corresponds to the case where the other singularity is at - oo. 
Figure 2-10a shows how to obtain the right limit point, U.- Draw- 
ing IN parallel to QO, because Z is now at O and take, arc Nee >. 
Drawing IU. perpendicular to Os, then results in U_- 


Figure 2-10b shows how to obtain the left limit point, U 


1 
Draw IM, parallel to QP, which is horizontal since P is at -o0o, and 
take arc M, 1 ). Drawing IU, perpendicular to OJ,» then results 


in U,- The geometric locus of U is the part of (U) circle between 


U, and UL: 


II-13 - Design Techniques 
In the preceeding section it has been derived a method for 
finding P and Z, given the desired location of point U on its locus, 

i.e. the magnitude and phase of QU, the unity vector. This method 

is restated below in a step by step form, Figure 2-8a. 

Step 1 - Considering Q as if it already were a point on the root 
locus, draw the vector diagram at Q for the uncompensated 
system and obtain QI, the "uncompensated" unity vector. 

Step 2 - Draw the circle diagram, composed of the (M) circle, of 
radius R= 1 and whose uppermost point is I, and the (U) 


2d 
circle centered at I and of radius r = Ao sin d. Fix the 
d 


oY 


limits of the locus of U on the (U) circle. 

Step 3 - Given QU as the desired unity vector for the compensated 
system, draw OJ penperpendicular to IU which cuts (M) 
circle at J. Draw angles JON = b. Then Z and P are deter- 
mined by drawing QZ parallel to IN and QP parallel to IM. 

In this section, it will be seen how the "U"' vector (i.e. the 
location of U on its locus) is selected to satisfy a particular 
condition. 

(a) Design for Minimum Root-Sensitivity. 

Since the U-vector is the scale used to measure the individual 
sensitivity vectors, the larger the scale, the smaller the magnitude 
of the sensitivity measures. Thus one possibility is to design the 
compensation for maximum magnitude of the U vector, that is minimum 
root sensitivity to open loop singularities. Figure 2-8a shows 
that maximum magnitude of QU is obtained by placing point U near the 
lower most part of the (U) circle or more exactly on the extension 
of QL. If such point is not within the locus of U, then it can not 
be a location for U and one must select the lowest point which is on 
the locus. In Figure 2-9a, this is point UL: Thus QU. is the selected 
unity vector, and with this given, one can proceed to the 3 step pro- 
cedure outlined in the beginning of this section. With such a design 


5.3 and 53) are all minimum, for all i's. 
(b) Design for Constant Damping When K Varies. 

Another practical problem is to compensate a system in such a 
way that when gain K varies about its nominal value, the dynamic 
response of the system does not change. This calls for a constant 5, 


i.e., a root locus that remains tangent to the radial line OQ at the 
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neighborhood of Q, Figure 2-8a. How can Z and P be found to obtain 
such a root locus? It is now shown that this can be done by merely 
selecting point U so that the unity vector QU goes through the origin, 
O, of the S-plane. In other words choose U so that Q, O and U are 

in line. If the locus of U does not permit such a choice, this means 
it is not possible to obtain a constant 5 about Q for the given system. 
One can then choose the best solution available, by taking the "U" 
location that is closest to a straight line with QO. 


The above statement can be proven very simply if one recalls 


Equation (2-44) that S 
if a5 
—— C-— 
Oo J 
S., 
The specification here is to force $41 to have the same direction as 
OQ. Thus making = move on a radial line when K varies. This means 
s. 
s.° must be equal to the phase of QO which is also phase of a 
But from Equation (2-63) 


Ss. 


J = Ss. x ! pas - 
Sp = S j - 0 (2-64) 
re) K 


The phase of an must be zero, this means that the sensitivity vector 
P 
s. a) 
Sp? must be on the U-vector or conversely, the U vector must pass 


through Py at the origin of the S-plane. 
S 


A faster way to prove this is that we show that Te always lies 
S 


a must be radial, 


on the U vector. In order to er constant, S 


thus U vector must be radial. 


II-14 - Limiting Behavior and Special Cases, 
The magnitudes of gain sensitivities can cover the entire range 
of values from minus to plus infinity. Yet, intuitive notion of 


“sensitivity'’ as a general concept in closed loop systems makes part 
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of this range unreasonable. One part of the problem is a direct 
consequence of the sensitivity definition, while another is associated 
with its first order approximation nature. A better understanding of 
both facets can be gained by an examination of limiting cases. 
In general, closed-loop poles depart from open-loop poles for 

low values of gain, and proceed to either open loop zeros or unbounded 
values as the open loop gain becomes very large. The gain sensitivity 
as given by Equation (2-13) could be written in the form 
s.J = 1 (2-65) 

n 1 mrtn 1 

a 2 as i = ees 
As K approaches zero, the closed loop root ee approaches the open 
loop pole from which it derives, i.e. Ste . Then the term, the 


al 


above relation reduces to 


"6 (2-66) 


Similarly as K becomes very large, n of the closed-loop poles approach 
open loop zeros. If the jth closed-loop pole is one of these, and 


it approaches the ith open loop zero so that 


Then, om 


J a 
Sy ep al ———— (2-67) 
K => co ( 1 ) 
2; - S5 
Finally, m of the closed-loop poles have no zeros to go to, and hence 


become very large relative to P, and Za. The sensitivity for these 


1 
poles is 
Ss, 1 3 
s J __——_—_—_— nm 1 -s j (2-68) 
& > ZiaP, 2 Ay ™m 
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Here m = (order of denominator) - (order of numerator). When the gain 
is sufficiently large for the open-loop zero db line to intersect 

the high frequency asymptote on Bode plot, the open loop transfer 
function is approximately 


L(s) = 7 


where m= (order of denominator) - (order of numerator). Then for 


the points on root locus 


1 + L(s) = 0 
Ss = -S, 
j 
or l + K = QO 
m 
( 54) 
or aT. a af EK (2-69) 


Thus, the sensitivity of the unbounded pole from (2-68) will be 
Ss, 7 
# >> Z,>P, m 
Equation (2-70) indicates that sensitivity increases as the m root 
of K, as gain is increased, although for the finite gains the sens- 


itivity is always finite. Referring to Figure 2-11, the above special 


cases can be simplfied as: s. 
J 

Ss, 2 0 
when K -> co K 

For poles going towards "0.L. zeros" cee 
when K ~~ QO K 

s, 
For poles going towards infinity Sve = j in general 
m 


Soi - "V -K when L = K_ 


m m 
S 


Another circumstance in which the sensitivity can become very 


large is revealed by Equation (2-45), and that happens when Q, (the 
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Fig. 2-11: Poles going towards open 
loop zeros. 





Fig. 2-12: Poles going towards infinity. 


au 


residue evaluated at pole TY) is not finite. This is to be expected 
since the sensitivity factors defined thus far have not considered 
multiple-order, closed-loop roots. As far as the gain is finite, 

an infinite gain sensitivity always indicates multiple-order closed- 
loop poles. 

A special situation of considerable interest can occur when a 
closed loop root lies between an open-loop pole and zero which are 
much closer to each other than to all other open-loop poles and zeros. 
This is the so called dipole case as in Figure 2-12. The sensiti- 


vity for the bounded closed-loop pole will be, approximately 








ate A ; wm (2478 4)(85-Py) 
K 
Z ; 8 it 8 : P “4 , 4 
i j j i 
(2-71) 
Differentiating equation (2-71) to find values of * which make, 
3 
sensitivity maximum, reveals that maximum value of s,? will occur 
8 
when = = bait Pi for which s_3 becomes 
2 
8 
J « 3 " 
Sy t (2, = P,) (2-72) 
max 


IIl-15 e- Sensitivity at Irregular Points 
Expanding the gain K = K(s) in series to include the higher 


order terms we have 


Ko = K +k Ae +1 d7K ps” + noe (2-73) 
° ds 2 ds- 


and then assume the first n-l derivatives vanish. This means then, 


AK = K-K = _1 dkAs" 
n: n 
ds 
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or 





As 8s 
8 /AK (2-74) 
i 
Since As is finite it is suggested to be substituted for Sy 


OK 


at irregular points. 

For any AK> 0, ds has n values separated by_20 which represents 
the outgoing branches in Figure 2-13. IfkK <0, n erties with the 
same separation as before are again obtained, but they are displaced 
from the previous set by Ti/n, that is the two set intersect. This 
could be stated as: The incoming branches, which at their junction 
represent an nth order closed-loop pole, are evenly spaced and sep- 


arated from each other by 27//n. The outgoing branches are also 


separated from each other by 27//n and are midway between the incoming 


branches. 


II-16 - Sensitivity Functions for Alternate Transfer Function Forms 
The transfer function could be written either in root locus or 


Bode form. For Figure 2-6, open-loop transfer function can be written 





as n n-l 
— C(s) ag s + as + e-e- + a 
E(s) Pe + b parael + eoe- b 
l men 
> at Tt 
a,s ‘ (s + Z,) 
K i=(0 =Xk iw] (2-75a) 
> b men-i ss p 
1° Tl. (s a 8 > 
i=0 i=] 
- 8 
Tr [e+ l 
= K' im] \“4 (2-75b) 
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The Equations (2-75a) and (2-75b) represent the root locus and Bode 
form respectively. 
The equations for the gain sensitivity are the same whether open- 


loop transfer function is in root locus or Bode fom, i.e. 


js j 4 
Se Po (2-76) 


This is because in comparing (2-75a) and (2-75b) it is obvious that 





TT 
‘ os ww 
K Ki Zs K x constant 
TT 
tuna 
and then, 8 ds ds 
ee 
K K' 


It is however necessary to modify the open loop pole and zero sensi- 
tivities for terms which are written in Bode form. For poles and 


zeros which appear in open loop transfer function in Bode form, i.e. 


S or S , the sensitivities are 
> 1 7 \+ 1 
i at 8 
i 
S 


j 
i I (2-77) 
i Z (2; 85) 
and s, = a 
Sp = j K (2-78) 
i Pi (s, ° P,) 


Frequently open-loop zeros and poles will occur as complex 
conjugate pairs, and variations in the system will change both zeros 


or poles. For example, consider a complex pair of zeros, Z, and Z 


i De 
which are defined by their frequency, w, and damping ratio S$, i.e. 


Z. eat SQ + JOVI = oe 


} 


(2-79) 
z. = 0 - jo¥l mG 


2 
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For this situation it is useful to define frequency and damping 


ratio sensitivities as 


soi = me Se + i 2 
Zz, 2d wz, 20d (2-80) 

oe ree 

S 3 62, =d§ Zz, aS 
Where ead and O84 are given according to (2-42) or (2-43) (In case of 

3Z1 320 dz, 2, 42, dz, 


poles). Components as To” ae? qb” and 7¢ can easily be eval- 
uated from (2-79). Then when the complex pair of zeros appear in open- 
loop transfer functions in root locus form, i.e. (8742 Sant 2), the fre- 


quency and damping ratio sensitivities for a complex pair of poles or 


8 + 5 Sj s q20s, Sej 
zeros are S19? < 20 DD <5 > iigt de + j_ > ee 
8, - 2508, + 8, - 23s, + 


Where the upper sign is to be used for zeros and the lower one for poles. 


288 the sen- 
wo * | 


If the term is written in Bode form, i.e. 8 
s|We 


sitivities are 5 3 48) (50 ai 8 4)8x 
w 
a 254s +0 
™ j j 
8 = J 
Sq" aeege ZEEE (2-82) 
2 2 
- 2 Ws dis u) 
5 ey 
For some cases it may be more convenient to define a complex pair in 
terms of their real and imaginary parts, i.e. 


Z4 = a+ jb Zo 2 a- jb 


In this case with the terms in root locus form rene 2as + (aoe: b*)), 


the sensitivities are 





a 
a 3a a2, da Z, da 
g 8 
. ot $3 
eee Clee 1) 
zy - aH Zo - os 
8 
= s.J 1 , 1 
(a+ jb - sj) (a - jb - 8) 
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54 2(a - 8,) as 
2 2 2 
a o oo a +b 


and in a more general form as in case of Equation (2-81) it could be 


written as s, 
* . = - ae 
oy = 2a8 a a +b (2-83) 
ee oe a 2 
a 2 <ep- a +b 


where the upper sign is to be used for geros and the lower one 
for poles. If the term is in Bode form (197/¢2 + b7)] + 


2as/(aé + b*)] + 1}, the sensitivities are 
{ 


s e 

nS ao Paa2. he Ps 

a > 2 , ; 2 s 
a +b s, = ae a +b 


(2-84) 
8 
2a - 8.) s7 


S 
j eg ( 
~ ” 2.2 ) a) 
a +b oy a seiae + b 


II-17 - Example 


In Figure 2-14, suppose the gain K, is given as 
s+ 
8 


K = - 


Find radius of curvature, and pole sensitivity. 


We can write aK. el ee 
ds 2 
$ 
are Fe 
ds* ~ 
ak _ 6 
de? ty 


Double points occur at dK _ 9 oF S=+ 1. Substituting S$ = 1 eb 
ds 
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Fig. 2-13: Outgoing Branches (n = 2) 





S=X+ jY 
Z2 ut jv 


Fig. 2-14: Example of 
Section II 
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we get 
dK 


& 
N 
) 
pete 
He 
SS. 
® 
Racal 
a 
“~&. 
ae 
e 


On circular portion of the root locus, 


dK 
4 Im ds 
2 


dx dk 
Re ae 


The pole sensitivity is 


8 2 2 
j = ds =] o es ele § 
Sx K OK ‘Willies 2 
les 
2 gs +i 
8 2 
i 
Sea s+ 8 
1 
BS 2 - 
8 


On the circular portion 


53 - 9 0, x. ‘ je bast J 


e€ - € 


The pole sensitivity at s = -l1 (a double point) is infinite and hence 


the suggested substitute, i.e. As is evaluated. Here n = 2 and 


VOR 


then 


As 


YY Ok 


e +1 evaluated at s = -] 





II! 
SENSITIVITY AND LARGE PARAMETER VARIATIONS 
III-l1 - Discussion 

This section deals with the most practical aspect of the sen- 
sitivity and is quite utilized in design consideration of the feed- 
back control systems. It discusses sensitity in two particular cases 
one with leakage transmission, and the other case when leakage trans- 
mission is zero. It relates the sensitivity of the system to return 
difference and null return difference. 

In previous sections it has been discussed that in most cases 
we are facing the variation in the values of an element in an 
engineering system. This element could be a vacuum tube, transistor 
or sometimes could just be a passive element. It is possible by 
means of feedback around the troublesome element, to achieve the 
desired reduction in sensitivity. Let the troublesome parameter be 
k, as in Appendix 1. We relate k to the controlled source S and 


control variable C, i.e. S = kC and, 


kt ,t 


ci os 
T = toi? ro kt. (3-1) 


with the fundamental signal flow graph as shown in Appendix 1. 


We defined the classical or system sensitivity as, 


Near Say , 
Si » on ke Ik /k (3-2) 


so sy is the relative change in T divided by the relative change in 


k for infinitesimally small changes in T and k only. 
Case 1: 


Consider the case where tis T, i.e. leakage transmission is 


j 
negligible. Tiis is a sitvation in many feedback systems. By straight 
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forward differentiation of Equation (3-1) we get 
gt ‘a QrT - sT/T -. k \T 
k aE dk/k 7 Ok 


= k 
ke. .£ eines 
ci os 
l - kt 
c 


(l = kt.) + tg (kt. t.) 


2 
8 G1 - Btast 


k(l - kt.) tae 
ee gee (l= kt y? 
cs 


and from that 


T 1 1 
aes “< - fort ,=0 (3-3) 
cs k 


where Fy is the return difference for reference k. It is clear now 


why it is insisted that all t.,‘s in the fundamental feedback equation 


ij 


should be independent of k. To have a small sensitivity then |F 


k 
should be large, i.e. (-kt. .) which is usually called the loop gain 








should be a large number. 

Since Cig and possibly k is a function of frequency, so is he 
and one must consider the frequency range over which a small sr is 
desired. It is interesting to note that as the magnitude of the loop 
transmission determines the sensitivity, its phase also determines 
the stability of the system. 

Case 2: 


When the leakage transmission to ¢ 0, then by straight forward 


differentiation as in Case 1 we get 


aT Be © oi os 
“ie a - ke)? 
cs 
. a Se 
T(1 = kt) 
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Fig. 3-1: Signal Flow Graph for Example (1). 





Fig. 3-2: Signal Flow Graph for Example (2). 


Input 1 a ra a 1 output 


Fig. 3-3: Signal Flow Graph for Example (3). 
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This could be written in the form 





st > i ta, as 
k (1 - kts) T(1 © kt.) 
and since from Equation (3-1) 
kt .t 
oi l - kt 
cs 
Pee, 
then T 1 oi 
Sk SD ke 1 i 
Cc 
t 
ne 1 : of (3-4) 
Fy T 


sy can also be written in terms of returm difference and null return 


difference. From Appendix 1, Equation (A3-3), we have 


Cp Ne Fe k 
g 
T FL 
then s" = -.— ° Sa (3-5) 
k k 


where Fy is the return difference with reference to k and Fy is the 
null return difference also with reference to k. Equation (3-4) 
suggests a practical mean to measure the sensitivity. It is noted 


that in case t = 0, Equation (3-4) reduces to Equation (3-2). The 


i 
applications of Equation (3-4) are illustrated by a number of simple 
examples. 
Example - 1 

In the single-loop system Figure 3-1, if the sensitivity desired 
is with respect to the forward transmission, the leakage transmission 
is zero, because with K = 0 there is no transmission from input to 


output. The sensitivity in this case is the reciprocal of the return 


difference, i.e. 


a3 


Sos = = (3-6) 


Example - 2 


In the same single-loop system, if the sensitivity of interest 
is with respect to the feedback transmission, B » the leakage trans- 


mission will then be K, and Equation (3-5) yields 


SB TE ( - pata): ar (3-7) 
In this example it is noted that when $= 0, i.e., there is no feed- 
back around k, then the leakage transmission toi =lxKxil1» kK. 
Comparison of Equations (3-6) and (3-7) demonstrates the well known 
fact that a value of K large with respect to unity results insensitive 
to changes in K, but with the closed-loop gain essentially equal to 
1 (when K is very large) the sensitivity of T with respect to B 


B 
approximates the value of -l1 = 1 /18° ? 


Example - 3 

Figures 3-2 and 3-3 represent alternative designs for a feed- 
back amplifier. Three amplifier stages are available, each with a 
gain a, and the overall gain to be realized is specified (less than 
a), Negative feedback is to be introduced to reduce the sensitivity 
of the overall gain to variations in supply voltage. To a first 
approximation this voltage variation can be considered equivalent to 
a variation in each value of a. The question to be answered is which 


configuration gives a lower value of sensitivity of the overall gain 


with respect to a? For Figure 3-2, the overall gain T, is 


1 
* 3 
T, *® a 
1 1 - £8 
The sensitivity of tT with respect to a is three times the sensitivity 
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of an individual stage with feedback, or 


1 3 
Ee Baa (3-8) 


For the second configuration, with a single overall feedback path, 


the overall gain T. is 3 
2 a 
1 - Pa 
wf 
and S 2 ~~ 3 7 (3-9) 
a 3 
l - Poa 


The two sensitivites (3-8) and (3-9) are to be compared on the basis 
of the equality of T, and 1. This equality establishes the relation 
ji 2s)? =. r- fa” (3-10) 


substituting (3-10) into (3-9) we get 


2 3 
a 3 
(i - Aa) 
Thus, the single overall feedback path results in a system with less 


sensitivity, to changes in a. 


III-2 - Large Parameter Changes. 

A serious shortcoming of Equations (3-1), (3-2), and (3-3) is 
that they apply only for infinitisimally small changes of k. There 
is therefore uncertainty as to their applicability for moderate or 
large changes in k. It is therefore found useful to use a new 
definition, defined as follows: 

Let ar 6 represent the nominal or original design value of 
the system transfer function and of the element under consideration 
respectively, and let Te Ke be the corresponding (final) values 


at the new value of k. Thus, 


af 


T. = Pract AT ke « kt Ak (3-11) 


The new sensitivity function is defined as 











t D 4t/T. lil 
k Ak/k . 
Case 1; 
When it is assumed that tui = 0, then from Equation (3-1), 
ores ee 
7, 7° Toke ag 
o cs 
k.t .t 
T «= _—f£ cios (3-14) 
2 _ Beck T.~ -T 
{ ° 
7 AT/T T 
sce Acie lnkr sok (3-15) 
f re) 
ke 
Substituting (3-13) and (3-14) into (3-15) we get 
l- T l - ka 7 k tig) 
gi - ae - k (1 r Kota? 
k kK. “ k, ke “ k, 
ke kK. 
l-ket F 1+ L 
o cs kA ° 


where F. is the initial return difference and L, is also the initial 
oO 
loop transmission. It is thus shown that when the leakage transmis- 


sion com is zero, the classical and new definition have the same 
value, although they are defined differently. With the new definition 
there is no uncertainty as to the effect of large parameter variation. 


Equation (3-16) is exact only if c is zero, but it may be used 


i 


whenever t. << T T° Fortunately ti is zero for a very large class 


i 
of feedback amplifiers. From Equations (3-11) and (3-16) we have 


T =f 
f Oo . fae Ak 
T 
£ Fa emma ace 
ln 1 Ke na cA 
le i ke 
( ‘. ) 
a 1_oco— 
l =k c ky 
Therefore i k : 
2-1-0-9(—7 
t o cs 
k 
r(temny a” ) 
£ ue o cs 
k 
and _o 
T) ke 1 ~ 
 * 1+ L (3-17) 
f re) 


where Lo a ek (tu is the original or the nominal loop transmission 
and Equation (3-17) relates the original and final values of the 
closed loop transfer function which is of very practical design 
application. 

Case 2: 


In previous cases we assumed t. was negligibly small in com- 


i 


parison with iS and T. over the range of variation of k. This is not 


f£ 
always true. Suppose in a system whose signal flow graph represent- 
ation is that shown in Figure 3-4 parameter D varies substantially, 
and feedback is to be used to reduce the system sensitivity to the 
variations in D. If we choose D = k of the fundamental feedback 


equation, then t. may not be small. However the problem can be 


i 
handled by method in Case 1, if we simply let k represent the entire 


element, i.e. k =  A(B + ED) 
1 - AED 
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Fig. 3-4: A Signal Flow Graph with Leakage Transmission. 





Fig. 3-5: Basic Feedback System. 





Fig. 3-6: Basic Unity Feedback System. 
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The element k is then incorporated into the overall system feedback 
configuration. The problem then reduces to Case l. 

In the above manner the majority of feedback problems can be put 
into zero leakage category. For those situations where a parallel 
path from the overall system input node to the output node (from I 


to 0) exists, the following results are obtained: 
Since T #& tot + KE stoa/( ° kts) 


using Equation (3-11) and notation (3-12) we get 


T AT/T. 1 (. a Ti 





S SB 2 


k Ak/k, (le k tag) 


t 
T = i ag _ of | > 
or Ss. ( 1+ L ) ( T, (3-18) 


Equation (3-18) can also be written as 


g(t oy) = 1 (3-19) 


1+ tL 
o 


and from (3-19) by using Equation (3-17) we get 








Ty o ta _ (k /Ke ) + L, 
Te & toi 1 + L 
or Ty : (1 + L) 2 (t 9, /T,) 0 © k /ke) 
T. L + k /kf 
o Oo ° 
ie, aye, MO ee) OL ke i 
T) Le i (ko /k e) 


III-3 - Further Investigation into the Meaning of Sensitivity Function, 


In prev ious discussions we defined the system sensitivity as 


T 1 
Sy = rs (3-21) 
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where L is defined as the open loop transfer function, and T is the 
closed loop transfer function. Let us now investigate the value 


dT/T 
dL/L 


Referring to Figure 3-5 we have 


T so a _ © L —___GH___ 
1 + GH H 1+ GH 
But L = GH and dL = HdG 
dT/T = _L dT » L dT 
dL/L T dL HT dH 
® 221 eS 
1 + GH 1+tL 


Now by definition dT/T D 56 i.e., sensitivity of the closed-loop 
dL/L Q 
to open-loop perturbations, and comparing Equations (3-21) and (3-22) 


we get uy 2 ¢C : 
Sr So (3-23) 
Also from relation 
T C 1 
Se) “ue rirean 


we can think of the sensitivity functions as the error transfer func- 


tion for a unity feedback system with forward elements GH, because 





R-C = £ 

or RGH 
Ro" CH s 

then E 7 eh , a 
R 1 + GH 1+ L 


In control system, the error transfer function is always desired 
small, ideally zero, and we may therefore observe that a small sensi- 
tivity function is desirable. 

Since the sensitivity function is dependent upon frequency and 


is a complex number for any arbitrary frequency we must further inquire 
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into its meaning. 


Let L rel? 
then s é. sr *, i 
1 + red 
= l 


(1 + Peos 9) + jC Pf sin d) 
or cc ae es 1+Tcos § ; vets d 


j 
> 1+ 2Icos fp +P? 1 + alcos $ +17 
(3-24) 
or 
5° Sx —— or -jtan"'] [sin § 
V1 + 2cos gd+f ~ 1 +I'cos @ 
(3=25) 


It is desirable to interpret both of the above equations from 
a physical view point. We then proceed in this direction. 
A - A Synthetic Interpretation of s 
a ee eS ee 


We define the following four quantities: 


gM D diT|/\|T — er unit increment in lr (3-26) 
M dq| LY /|L| per unit increment in | L]} 


d/T | 

s* 4 ——— = {increment in fot (3-27) 
d / tL increment in / L 

M D diTI/ITI _ _per unit increment in |T | 

>i. d/u increment in / L Se 

st z d [T = increment in mt (3-29) 
ai LI/\L\ per unit increment in |L| 


Where |a| and /a indicate the magnitude and phase angle of the 
complex number a, respectively. It is necessary to compute |T | 


and /T_ in terms of | L| and /L_ in order to evaluate these functions 


L ® Pr _iP 


Therefore |u| e |! and /L_ #® g (3-30) 
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1 (39 


and T = a T+ Tie (3-31) 
let H = | nj eI and G = |cle/e (3-32) 
therefore 





(* al ¢ - tan”? fring fring - 18 
lH fy feos 6 +22 ~ 


(3-33) 


T = 


We are now in a position to compute the functions defined by Equations 
(3-26) to (3-29). Note that 
[ = | Hl @ Ie a" «= | aledic| (3-34) 
dQ = p+/o df = d fe (3-35) 


From Equation (3-26) we have 


x20 di Ti /\it\ 
M dILi/IL \ 


= PHI AD y 2 cos g +f 


But |H| is not a function of the variable part of r then 


Ss = V1 + 20% cos § +17 | 1+ 2! cos +f? 
df 


1 + C cos g 





1+ 2%cos g+!? (3-37) 
By similar technique we find 
* eee cy (3-38) 
™ — ae (3-39) 
s 7 ——-[sing (3-40) 


1+ 2%lcos g+f[? 


From above results it is shown that: 
¢C et 
Re Ea = Sy S 
¥ 


=| 


So gM sf Sy? and ee are all real mumerics which are desired small, 


ideally zero. To imvestigate the utility of the sensitivity function, 


(3-41) 


eK Sw 


(3-42) 


the mapping of the loci of constant s and si on the polar plane 


(hence all four) will be considered. 


B i 1 ai o£ st beet’ onto ais peas goes 


4h 





We require 
1+ [cos 6 


= M (3-43) 
142%cos f+! 


Where M is a constant. Changing to cartesian coordinates we get 


Es oo (3-44) 
1+ 2xee+ x + y 
or 2 
1 2 2 h 

x+ (1 - a } + y ( 2] (3-45) 
Thus the loci of constant S. are a family of circles with center 
coordinates ve: og - [wate y 20 (3-46) 
and radius = a (3-47) 
: pe-Bepying apes ES Eee onto the sei rian 





We require 


___clsing ssi, Q (3-48) 


1+ 2fcos ¢ +f? 


where Q is a constant. Changing to cartesian coordinates we get 


Q = Sy (3-49) 
l+ 2x+x + y 
oF (+1)? +(y+1_.)* = /1_\’ (3-50) 
20 20 
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Thus the loci of constant sh are a family of circles with center 


M 
coordinates x = el, ys e-il (3-51) 
2Q 
and radius = 1 (3-52) 
2Q 
Both families of circles derived for sy and si are plotted in Figure 


3-7 and it is observed that the region of poorest sensitivity is 
near the (-1, 0) point, the sensitivity improves in every radial 
direction from that point. 

D - Mapping of s. Loci onto the Polar Plane. 


From Equation (3-25) we require 
1 


Vi+2%eos p+? 


where M is a constant. In cartesian coordinate it is 


(x41)? +y? = (+) : (3-54) 
M 





= M (3-53) 


The loci of constant | Ss. | are a family of circles with center 


coordinates, x = -l, y = 0 (3-55) 
and radius = ji 
M 


c 
E - Mapping of [o Loci onto the Polar Plane. 





From Equation (3-25) we require that 


ol ( sin , 
tan i+ 1 cos p = C (3-57) 


or in cartesian coordinate 

y = tanc' (1l+x) = C(1 + x) (3-58) 
The loci of constant phase angle of the sensitivity function are 
therefore straight lines through the point (-l1, 0) with slope C. 
This family of lines comprises the orthogonal trajectories to the loci 


of constant | s 


{[1Jey 


F - Result of the Interpretation. 


For any point on the polar plane 
2 





(3-59) 





where each term is evaluated at the point in question. Thus the 


magnitude sensitivity function is a vector measure of Sy and si at 


every point of the polar plane. 
The only extent to which the phase angle loci will be o¢ use 


comes from the consideration of the phase angles zero and™ rad. 
2 


Along each of these lines the imaginary and real parts of so 


tively are zero, We conclude that along the phase angle locus corr- 


respec- 


esponding to zero rad., magnitude changes in L cause only phase changes 
in T. The reverse is true along the phase angle 71 rad. locus. 

Since frequency response design work is generally carried out on 
the amplitude-phase plane, the mapping of the magnitude Fensttavies 


function loci onto this space will] now be considered. 


III-4 - Mapping of the Magnitude Sensitivity Function Loci onto the 
Amplitude-Phase Plane. 


The magnitude sensitivity function loci can be mapped onto the 
amplitude-phase plane simply by inserting the variable change 
Z «= Ly | 


into Equation (3-53). There results from this substitution, 


ee? 4 20° cue QO+ 1. (2) 2 (3-61) 
Plotting this on the Z-9 plane gives the ide via, 
The ultimate purpose for the creation of this mapping is to use 
it as an overlay in design, just as is used in the Nichols Chart. By 
a simple manipulation it will be shown that the Nichols Chart itself 


is already calibrated for our use. In previous sections we concluded 
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Is} 


queysuoD 


ue) 


so So = a (3-62) 
where, L(s) = G(s) H(s), Figure 3-5. It is intended to point out 
that this sensitivity function is obtainable directly from a convent- 
ional Nichols Chart, given the open-loop frequency locus 

L(j) = G(j®) H(40) 
for the linear system in question. The procedures is quite direct. 


Both numerator and denominator of Equation (3-62) may be divided by 


L(s) gi : ing b 


1+ (8) 
We seek the polar coordinate mapping of lines of constant magnitude 


sensitivity function given by 


or 
Le) = M (3-64) 
1 + L(s) 
let L(e) | «© @+ fy (3-65) 
therefore 
frst re ee (350m) 








and from Equation (3-66) we get 


uz =—«\ 2 2 M : 
M -1l M -1 


We recall the equation for the M circles of conventional design 
(ref. 9) and note that it is identical to Equation (3-67). Thus we 
have the interesting result that the M-circles are to the open-loop 
frequency response curves as the magnitude sensitivity function loci 
are to the inverse open-loop frequency response. The Nichols Chart 


is just a plot of the M-circles (as well as the N-circles of constant 
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closed-loop phase) on the amplitude-phase plane, therefore the mapping 
of the magnitude sensitivity function loci onto the maplitude-phase 
plane is the Nichols Chart, provided that we read intersections as 
derived from inverse open-loop frequency response curves. 

The combination of the logarithmic amplitude scale and linear 
phase scale on the amplitude phase plane, with the symmetric Nichols 
chart is shown on Figure 3-9. Consider a typical open-loop locus and 
its inverse, on the extended Nichols chart as in Figure 3-9. The 
inverse function ja) is obtained on a@ point by point basis graphe- 
ically and is symmetric with the direct function L(j)) with respect 
to the point (O DB, o deg). The Nichols chart is shown over a 540 
degree range and is seen to be symmetric (Ref. 9) about the axis 
phase-0 deg. 

One further simplification may be made as shown in Figure 3-10 
If we imagine that the phase calibration for u (40) locus on the 
Nichols Chart is opposite to that commonly used then the same Me- 
contrours used in ordinary control system design, measure sensitivity 
function magnitude for use to the extent of being properly labelled. 
In this case L( ji) and Lt ju) are symmetric with respect to zero DB 
line, with phase measured in opposite directions for each. For the 
particular example of Figure 3-10 we see that the maximum magnitude 
sensitivity function value is M = 2.0. 

If it was desired, the magnitude sensitivity function for the 
closed loop system as a function of frequency could be read directly 
off the Nichols Chart as the intergections of the M-contours with 


the inverse open-loop frequency response curve. 


a 


The simplification introduced by foregoing interpretation of 
the sensitivity function magnitude and phase loci could be very use- 
ful in the design of linear systems with regards to sensitivity 
specifications. 

Special Case: 

If the forward elements can be separated into two tandem parts, 
one of which displays sensitivity to the environment (denoted by an 
"environmental parameter" x), viz., 

G(s, x) = G, (s) Go (s, x) (3-68) 
Then we may compute the sensitivity of the closed loop locus of the 


overall system to the environmental parameter (denoted s") according 





to 
sf . dT/T a dT/T dx/x 
0 dL/L dx/x  dL/L 
i x 
= S$) 8% (3-69) 
14 G& 1 CG 20 
or s = 8 x Ss, 8 (3-70) 
oO 
In the present example 
so aL/L | dG, (8,x)/G,(8,x) 
x dx /x dx/x 
Co 
i S.. (3-71) 
Therefore G 
T C 2 
Ss. * 55 %& (3-72) 


where the first term of the product is evaluated graphically, and 

the second term is found by differentiation. It should be noted that 

if the parameter x appears in G,(s) as G, (8) = xG,(8), then 

Be. = 1 and 5° ze 5° 
x 0’ 


T 
So = Sy » where k appears as in Figure Al-l of Appendix l. 


which is in accordance with our relation 
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Fig. 3-9: Typical L and its Inverse on Extended Nichols Chart. 
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Fig. 3-10: L and he on the Nichols Chart. 
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III-5 + Philosophy of the Frequency Response Approach to Sensitivity 
Problem. 


Case 1 - When ¢t es QO, 
oi 


In previous sections we discussed how sensitivity and frequency 
are related and how by utilizing the Nichols Chart, the magnitude and 
phase of the sensitivity function and its peak magnitude can be obtained. 
Another approach which could be of high importance in most design cases 
is discussed as follows: 

Sensitivity specifications could be given in terms of the systems 
behavior on the jb axis. It is possible that the allowed variation 
in T(jQ) may be exceedingly small over one frequency range and fairly 
large in other frequency ranges. The specification may obviously be 


in a wide variety of forms. For example the nominal T,., may be as 


0’ 


shown in Figure 3-11, with bandwith ,, amplitude peaking Mo» and the 


1? 
specification may dictate that despite parameter variations, the 
magnitude peaking should never exceed 30 percent, the bandwidth should 
be in the range 0.86, to 1.20,, and that from zero frequency to 
0.50, \T| should never be less than 0.95. Finally, some thought 
should be given to the intermediate and the far off frequency ranges. 
For example, the behavior shown by dotted lines in Figure 3-11 is 
usually intolerable, as it results in high frequency ringing. A 
reasonable statement might be that for 1.20, KCOoOK< OQ, ‘ | T( 40) < 3DB 
and for W> W, T(j®) must not exceed a predetermined value. 

The basic philosophy of the frequency response approach to 


sensitivity reduction is deduced from Equation 3-17, rewritten here 





z. (k /k,) + L 
° o f ° 3-73) 
T * Teyw ( 
a © 
This equation applies whenever the leakage transmission wi is zero 
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and therefore could be used in most control problems. Suppose that 
at some specific frequency s = Ju plant parameter variations are 
such that (ko /k e) at s = ju, may lie anywhere inside the region 

indicated in Figure 3-12. Suppose also that - L,(j4) is given by 


the complex number A + jB located at Q. Then Equation (3-73) states 





that T oy 
e) - QV 

dt —> 

f QN 


The range of variation of the vector QV fixes the range of variation 

of T/T. at 5s = ja. Suppose for example, that it is required that: 
1.20 > T/T. > 0.80 

at 6 = jo. The problem is then to find the locus of “L542 which 

berely satisfies this inequality. A suitable procedure is as follows: 


Suppose we seek the locus of the boundary of Ho 


T 
of any fadius R (Figure 


> 0.80 with 








point N= 1+ ie as center, draw e circle c 


3-13). Next draw a circle C, of radius 0.5K, using as center the 


point (the boundary of ko/k which appears to be the closest to the 


Py 
first circle - this appears to be point J in Figure 3-13. The inter- 
section of the two circles determines a point Q. As a final check, 
use Q as center, and check whether a circle of radius 0.8R cuts ko /Ke 
locus only at J. If so Q is on the boundary of the locus of 

| T)/T,| 2 0-80 at s = jo. The above construction is repeated, using 
different values for radius R, until the locus of Q is obtained. Such 
@ locus for 1.20 > | a le 0.80 is shown in Figure 3-13. Clearly, 
“L Gia) must be located inside the indicated region in Figure 3-13. 
On the other hand if there is the additional requirement that the 


phase of T/T. must not change more than 20°, then the permissible 


range of L654) is reduced to the smaller region shown in Figure 3-13, 
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It is extremely important to note that due to Equation (3-17) the 
sensitivity requirements essentially determine the gain and band- 
width of L(s) with hardly any need to refer to the specific feed- 
back configuration that is used. At best one configuration may be 
more efficient than another in the practical realization of the 
required L(s). The latter is actually a problem in network synthesis. 


Case 2 - when toi ¥ 0. 


In case 1 we assumed that in the fundamental feedback equation 


tT - t.. + kt to—/@ © kt .)s the leakage transmission ty was 


oi i 


negligible in comparison with T, ower the range of variation of k. 
As previously discussed this is not always the case. When the leak- 


age transmission t, was not small in comparison with T, then Equations 


fi 
(3-18) and (3-20) were resulted which are rewritten here: 





" : 
T AT/ f a 1 ¥ oi ig 
ag Ak/k. 1+L, (1 Ey ) (3-74) 


TS =¢ (k/ke) + Ly 


T.-t - i¢ tL (3-75) 
£ © 


<4 (1 + Lo) 2 (t.,/T,) (1 - k /ke) 
T L, 4 k Ke 








(3-76) 


where Te is the value of T when k has the value k and 


4T = T,°T, 


£ 


Ak = ae as 


i is a transfer function which is independent of k. If tog is 


known and the desired T is known, then one may relate the maximum 


Now t 
re) 


permissible range of variation of T, over amy frequency range, to a 


corresponding maximum range of variation of T = ta: In other words 
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(a) 







“11 - ko/k.) 1 k /k, freq. locus 


tito freq. range 


(b) 


-to,/T, (1 - k/Ke) freq. locus 
= (es 
L 6 ) 

i KM 


(c) 


locus k ke 


A 
(oy ill 
1- t ,/T, (l= kjfke) 


Fig. 3-14: Construction for finding T/T.» when t , is not zero. 
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one can work graphically with Equation (3-75) in exactly the same 
manner as Equation (3-73) was used in Case 1. The only difference 
is that one must first decide what is the permitted range of vari- 
ations of T. - t ,, rather than that of T.. 
£ oi f 
The above is satisfactory if tot is known. On the other hand, 


there may be situations when t . is not precisely known until part 


ot 
of the design has been done. [In such a case some cut and try is 
inevitable. The following modification of the graphical procedure, 
suggested by Horowitz is an alternative method of handling the non- 
zero leakage transmission problem. Equations (3-75) to (3-77) are 
used, Sketch the frequency locus of tit, (Figure 3-14a) in the com- 
plex plane, over the frequency range in which control of variation 

in T is desired. Also sketch the ioci of k /Ke- From these two 
figures, obtain the loci of -(1 = k /Ke) (t/t) as in Figure 3-14b. 


Finally combine the latter and the sketch of k /Ke into one figure 


(Figure 3-14). In Figure 3-l4c we have 





in accordance with Equation (3-76). The design problem is to select 


the loop transmission L, 80 as to satisfy the specification of T/T: 


III-6 - Use of Polar Plot in Studying the Sensitivity. 


The relation, gt S © > 1 
K O 1 + L (5) 


suggests that the inverse polar plot could be useful in studying the 


sensitivity function. Above relation could be written as 


L (8) 
© Ll + L, (3s) 
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Referring to Figure 3-15 it is noticed that if OP represents L.'Go), 


then BP = 1 + Lo) and, 





ol 
L_  (j) 
S = -so = = = Med (3-79) 
1 + L, ( jW) 
As discussed in Equations (3-63) to (3-67) we have 
2 2 2 
x = 5 + y? iad —s (3-80) 
l-M 1-M 
where x = real part of L. (48), 


y = imaginary part of ge j®). 





Equation (3-80) represents a circle having radius = = 7? and 
1-M 


center at - 2 on the x axis. (3-81) 


The loci of constant M described above is useful for determining the 
value of M ax for a control system for which L, (8) and eventually 
L,(e) is known. The radius of the M circle and the distance from 
the origin to the center of a given M circle are both a function of 
M alone. Referring to Figure 3-16, where P is the point of tangency 
from the origin to the M circle, one observes that the ratio of BP 
to OB is a constant for a given M. 
BP = radius of Mcircle = m/or - 1), 


oB = M?/(1 - M°), 


Rae. Foss 
80 2 
sin s i - + 
M 
1 -M 


Then sin ‘Y is solely defined in terms of M. 


Oc OP cos Y OB gant ¥ 1 “2 - 1 
—_—- 3s i s lo —_-s oOo 
OB OB OB M2 m2 
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Complex Plane 





Fig. 3-15: Conventional Form of Complex Plane Diagram for 
Control System of Equation (3-78). 


Fig. 3-16: Geometry of M Circles. 


83 


or OB/OC = W/O = l). 
Hence C represents the location of -1 + j0, point to the same scale 
for which BP represents the length M/M” - l. 

For a given L.(s) function to have a maximum M value of a it 
must also be tangent to the M circle having a value of M . By using 
these facts the following procedure can be formulated. 

1 - Draw the locus of L”*(Jo) function on the complex plane 

with the value of K = 1, K being the gain. 


2 - With the value of M_ specified, determine the angle | from 


Y = “0 (at) 
att 


3 - Using the angle Y determined from 2, draw the line OD on the 
complex plane. 
4 - By cut-and-try, draw the circle that has its center on the 
negative real axis and is tangent to both L*(4e) locus 
and to the line OD at some point P as shown. 
5 - Draw a perpendicular PC from the point of tangency P to 
the negative real axis. 
0C(1/K) = 1.0 
then K = OC 
This value of K is the gain associated with L_(s) that will yield the 
specification of M.. Knowing the proper net gain K, one may then 
change by K the scale of the plot used in preceding steps and evaluate 


the sensitivity. 


IV 
SPECIAL TOPICS 
IV-l - Discussion. 
In previous sections we discussed the sensitivity function and 


showed that: 


Sr = a pase toe (4-1) 
k T 


where Fi. is the return difference with respect to the parameter k. 
The parameter k could be any parameter such as gain K, pole Pie or 
zero Zs and Co 
and output, as shown in Figure 4-2. Although the sensitivity function 


is the leakage or direct transmission between input 


has been defined and several of its properties examined, it is not yet 
apparent how one may utilize it in system design, or more fundament- 
ally on what basis its characteristics should be chosen. Similarly 
in order to compare the relative sensitivies of different parameters 
it is convenient to have an overall measure of sensitivity which is 


independent of either time or frequency. 


Iv-2 - Specification of the Sensitivity Function. 

A common design approach consists of making the “loop gain" or 
return difference, extremely large in the band of frequencies in 
which the input signals are expected to be. Then the open loop 
response is caused by appropriate shaping to fall off as rapidly 
as possible to unity gain outside this band. A maximum roll off 
rate of about 33 db per decade permits unconditional stability to be 
achieved, while a greater rate results in conditional stability. Thus 
realization of low response sensitivity by making the magnitude of 


the sensitivity function small in the desired band of frequencies 
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-1 


input A(s) B(s) M(s) 1 Output 
R(s) C(s) 


Fig. 4-1: A Basic Feedback System. 


- (3¢5) (<8) /K}) 
(a) 










A(s)B(s)(M(s)/K] 
R(s) 1 C(s) 


-B(s){M(s)/(S + 2,)] 
1+ B(s)(M(s)/(S + Z,)] 







(b) 


A(s)B(s){[M(s)/(S + Z3)] 
(1 - B(s)[M(s)/(S - zs) 





R(s) 


SA(s8)B(s)(M(s)/(S + Z,) 
I+ SB(s)(M(s)/(S + P,)] 


Fig. 4-2: Alternate System Configuration of Figure 4-1. 
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(c) -1 


S + B(s)(M(s) (S + P,)] 






P 
(S + B(s)] [M(s) (S + P?) 


R(s) C(s) 
A(s)B(s)[M(s) (S + P,)] 
S + B(s) (M(s) (S + P,)] 
(d) 
Input Outy: 
R(s) C(s 





Fig. 4-2: Continued. 
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inevitably results in the problem of system stability. In the 
following several criterions for the specification of sensitivity 
functions are discussed. 
(a) - The Forced Response Criterion. 

Suppose the desired response in the time domain is given by 
c(t) and the response error by e(t). Referring to Figure 4-l, if 


(8+ 2Z,) (8 + Z,) ---- 


M(s) = KG?) GFP) oe (4-2) 
then T(s) = C(s) A(s) B(s) M(s) (4-3) 


R(s) | 1+ B(s) M(s) 
and assuming variable element x is associated with M(s) only, then 


s'(s) = x qT(s) 4. x dM(s) /dx 


T(s) dx M(s) 1+ BCs) M(s) 
(4-4) 
for x = K, 2 or P, Equation (4-4) becomes, 
Y 1 
S(8) " TS BG) Mey i 
Zz 
T 1 
S09)" Te BG) May 8+ Z, ‘SS 
. Ne 
T i 1 
Sp(s) = SFP CTF BGs) Me) — 


The basic problem is determination of criteria for the specifi- 
cation of the sensitivity function for an arbitrary input. After 
the form of the desired transmission function has been ascertained, 
based on some dynamic response criterion, Equations (4-5) to (4-7) 
reveal that the system parameters which remain to be specified are 
the zeros of the sensitivity function. The poles of the sensitivity 
function and the transmission function are identical, except for the 


added poles Z, and Pi» provided the product A(s) B(s) M(s) in Equation 


j 


(4-3) is chosen so that no poles are added to those of the reciprocal 
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of F, Cs). The forced response criteria determines the zero of S(s) 
so that e(t) behaves in the desired manner over the interval or 
intervals of interest in the time domain. A considerable literature 
exists whereby one may relate time domain responses to the complex 
frecuency domain characteristics of a given network configuration. 
These results are useful when the inputs are restricted types and 
when the system poles are simple. Functions of third or higher order 
with multiple poles and possible poles on the j axis are difficult 
to deal with. If the sensitivity function is to be specified only on 
the basis of the forced response to a given input, then the zeros of 
S(s) are chosen as follows: 


Lf &(s) consists of a number of inputs of the for, 


R,(s) 2 ee e+e, (4-8) 
then the system output is given by, 

c.(s) = ‘are, MO (4-9) 
and the forced response is 

c(t) = a,T(-s,)e°j° (4-10) 

j j j 
Taking the derivative of © ,(t) with respect to a variable parameter 
x we get, 

ie. ad oe (4-11) 

dx j dx 

But from Equation (4-4) we have 

—ts = s° (-8,)/s (4-12) 
Then substituting Equation (4-12) into (4-11) we get 

—— : — T(-8,) S)(-8,) @ a A ac 
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Thus, if the sensitivity function s(s) has a zero at the pole of 
the transform of the input function, the system forced response is 
invariant with changes in x. If Jie are expected to lie in some 
known range of values, the sensitivity zeros may be distributed 
accordingly. 

The above result is also valid when 8 , is imaginary but in that 
case both the steady state amplitude and phase responses are of 


interest. For a sinosoidal input, the forced response and its deriv- 


atives with respect to x become 




















a 
c(t) ~ r T(js_) sin E + b cs) | (4-14) 
r 
where =| ImT( js _) 
#(js.) = tan Rete) 
and 
de (t) a | t( 5s ) 
r ig d r 
Sac? c (st + 9) =a 
ap(js_) 
T(js_)| —3>——_ cos (8 t + 0)\ (4-15) 
Writing T(js_) -|t038) JP (4-16) 





jo 
Then s"(js_) h, x dj Itle-~_) 


ee 


x dlt| dg 
- Tt 7 ax * J “ax 


= s.(js.) + §8,(j8_) (4-17) 


The R and I subscripts in the last equation refer to the real and 
imaginary parts of S(js_), respectively. Thus comparing Equations 
(4-15) and (4-17), if zeros of S(s) are assigned at + js_, then both 
S. and S, are zero and d¢_(t) becomes zero and is invariant with 


R I by 


changes in x. 
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In this section we domernstrated how variations in the forced 
response may be controlled by suitably specifying the sensitivity 
function. 4A design based only upon the forced response is likely to 
result in an unsuitable selection of the zeros of the sensitivity 
function, if the system is of high order and the input is of a known in- 
variant form. For example if a fifth order system is to be subjected 
to a sinosoidal input, then five sensitivity zeros must be assigned. 
Two zeros are given values corresponding to the poles of the input 
transform, there is no reasonable basis for selecting the remaining 
three. Furthermore, if interest is centered upon the part of the 
system response where the transient has an appreciable value, then 
even the forced response criteria outlined have little value in deter- 
mining suitable sensitivity zeros. Then another technique besed on 
“minimum mean square error” was proposed. 

(b) - The Minimum Mean Square Error Criterion 

The winimization of the mean square difference between a desired 
and an actual response is fully discussed in literature. This tech- 
nique is directly applied by Mazer to the sensitivity problems. Its 
main advantage is that the problem is solved entirely in the complex 
frequency domain, thus the necessity for obtaining the error in the 
time responses of high-order feedback systems in terms of the unknown 
zeros of the sensitivity function is avoided. Before proceeding with 


details, the following definitions are required: 


T(s) = —oe (4-18) 
N 
—~ oS ; 
S(s) D(s) (4-19) 
At(s) = 2- 1%s) s(s) (4-20) 


E(s) # R(s)A T(s) (4-21) 
Where R(s) is the laplace transform of the input and A notation is 
used to indicate a small finite, rather than differental, change. 


The mean square error is given by 


T 
as Lim — J e7(t)dt (4-22) 
T-> 00 or 


Here e(t) is the inverse transform of E(s). The square of the error 
is integrated over a time interval of length 2T, its average is 

taken and the time interval is then allowed to become infinite. It 
should be noted here that T(s) is the transfer function of the system, 


while T is the time interval. As shown in reference (2), 


joo 
e = 5 J” ArGs) AtC-8) §,(e)4s (4-23) 
ja 





Where p(s) is the power density spectrum of the input signal. Sub- 


stituting Equations (4-19) and (4-20) into (3-23) we get 


j@ 2 W(s)N.(-8) 
2 1 Ax ie) 0 
wt J 4 Deayp(-ay TESEC-#)P, (248 
(4-24) 





To minimize the mean square error ef one takes its partial deriva- 
tives with respect to each unknown coefficient of N, (8) and sets the 
results equal to zero. This yields a number of linear algebraic 
equations equal to the number of unknown coefficients. Furthermore 
if No(s) is of even or odd degree, all coefficients of odd or even 
powers of s must equal zero. In order to obtain a suitable density 
spectrum for use in Equation (4-24), aperiodic inputs such as steps 
and ramps, must be modified by being repeated periodically. Although 
the system is then designed on the basis of this fictitious input, 


the periodicity of the repeated input signal determines approximately 
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the interval of time over which the mean square error, for the case 
of an aperiodic input is minimized. 
Example: 
Design a feedback system in the configuratior of Figure 4-1, with 
the following closed loop transfer function 


5 
(o- + oF Th er 5) 





T(s) = 





(4-25) 


Since the system is of the third order, No(s) may immediately be 


written N.(s) = 3° +. a8 + a,s + a (4-26) 


If the variable <etement is the forward gain «, the mean square error 


becones: 
eS 2 
cea 1 y \ kK T(s)T(-s) 2 2 AM 2 2 
. 27 j ( £) Draji(-ef | 2” 2g eee ee 
= ja 
f,(s)4e] (4-27) 
Since a? is the sum of two integrals, one may solve for a, and a 


separately. Taking partial derivatives with respect to these 


quantities, 


lak - _T(s)TC-s) | 8)T 


Ky Sere 








» 42 
(4-28) 
Ene a) sf (e)ds = 0 


and similarly 


re 1 yi Ax\? T(s)T(-s) 
ea, MG ZG NK] CDG 8) 


Oo 
2 

| 2a, + a, 9, (0008 = 0 (4-29) 

Solving Equation (4-28) and (4-29) for a, and a, homogeneous set of 


equations results, the solution which is 


To solve for a; 
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- joo 2 
be 1 Qk T(s) T-s) 
da, = 2 j ) (- E.| D(s) D)-s) fe 2° (s “40,) 19, (s)ds e 0 


2jo 
(4-30) 
and from that 
1 He T(s)T(~s) 2 2 
273 Disyo(-s) (© (-8°)P, (e)ds 
a, = —<— oe (4-31) 


1 J a(8)T(<8)_ 
tr S ONES m 5 (3)(-8)9_(s)ds 


Thus a, is independent of AK, provided the latter is small enough for 
the first order approximation AT(s) given by Equation (4-20). Further- 


more p(s) being positive, a, is always positive and the integrands 


1 
are written as conjugate products in Equation (4-31) to emphasize this 
fact. The character of the input remains to be specified so that 

és) may be evaluated. Suppose the input is a sinusoid of frequency 
f = 1/2 rad/sec. and that the mean square error of the system response, 
over one cycle of the driving function is to be minimized for small 
changes in k. To generate an appropriate repetitive input signal let 
the input sinusoid be multiplied by a square wave of frequency 

- = 1/4 rad/sec. and unity height. The Fourier series and power 


spectrum of the square wave are 


1 1 3 l > 


wean? = - (sin (t+ 5 sin ft + 5 sin Zt + --- ) 
and 2 
Iq) - : (4) cc : 2) + §CO+ ZI 


+ 5{5(W 2) + (w+ 2)} 


* 213(= 2) + §(O+ 2) 


+ cone-e-ee (4-32) 


34 


Here 3( + n&) is the Dirac delta function. If this square wave 
is driven through a fictitious transfer function 


is) =. ———— (4-33) 


2 


then the power spectrum of the input is given by 


2 
P.O) = 9.) FG) 


G2 
7 Teo $4”) (4-34) 
for = (2) 

With #4) given in Equation (4-32), the intergrals in Equation (4-31) 
are evaluated by summing twice the values of the integrands at J = a 
2 ,? » etc., since they are even functions and thus always positive. 
The summations are continued until the desired accuracy is achieved. 
The values of the integrands fall rapidly with increasing 4) , since 
the degrees of. the denominators are much larger than those of the 


numerators., Using only the first five terms in the Fourier series of 


Equation (4-32) and substituting T(s}), D(s} and ¥ (s) we get 


J _T(s)T(es) _ 8)T(-s 3 3 =2 
i MESS | D(s)D(-s) (s) (-8) f (8)ds = 0.751(10 °) 


(4-35) 


i” _T(s)T(-s) yi v4 -2 
3 4 D(s)Dtcey 2) (8) @.(s)ds = 1.238(10 ") 
(4-36) 


” _T(s)T(es) | > 
Ee 2” Agee D(s)D(-s) (8)(-8)6 (s)ds wm 2.511010 *) 
(4-37) 
Then by Equation (4-31) we get 
a, = 0.404 


and the mean square error is 
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j 2 
2 (s)T(-s) I 2 
eo * ay 2 (4) Dayo(=ay Er# Ce Fa)" 1f,, Codes 


(4-38) 
Equation (4-38) can be written as 


2 i” 2 _1(s)T(-s) 
ot d. CE) D(ayDUcay (v8 CDP, (6)48 





| ee 2 
a - 2 





2” (4By Oe (-8)(s)a, “4 _(8)de (4-39) 


and utilizing relations (4-35) to (4-37) we get 
2 
e* = 0.00149 (48, (4-40) 


The sensitvity function is then 


N (s) 2 
T re) s(s_+ 0.494) 
= D(s) geaat s+1)(s+ 5) | 


The above techniques show how complicated the computation could 
be and in many cases the use of computer is the quickest way. The 
purpose of the above discussion is not a complete design technique, 
but to show how different approaches are made to the problem. 

For more insight into this problem reference 2 and 11 are 


recommended, 


IV-3 - Equivalent Input Perturliation Signal 
As discussed in Appendix 1, for any linear system, the transfer 
function T(s,x), relating a response to an input function can be 


expressed in the following bilinear form in terms of a real parameter 
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x of the system and the polynomials A(s), B(s), C(s), and D(s) 


T(s,x) = _P(s8,x)_ e _A(s) + xB(s)_ (4-42) 


Q(8,x) C(s) + xD(s) 
where P(s,x) = As) + xB(s) (4-43) 
Q(s,x) = C(s) + xD(s) (4-44) 


From Equations (4-43)and (4-44) we have 


A(s) = P(8,0) (4245) 
B(s) = —> [P(s,x) - P(s,0)} (4-46) 
C(s) = Q(s,0) (4-47) 
D(s) = — [Q(2,x) = Q(s,0)} (4-48) 


With the aid of these equations the logarithmic sensitivity function 
may be expressed directly in terms of the polynomials of the original 


transfer function. 


T OnT(s,x x 3T(s,x 

Tr ary wee Re em fg Ni 

But 
x - —* Q(s,x) ms 

T(s,x) P(s,x) (4-50) 
and 

dT(s,x) . > E + spied |) 

x Ox C(s) + xD(s) 


_B(s)C(s) - D(s)A(s)_ 
[a¢s,x)}° 


Substituting Equations (4-50) and (4-51) into (4-49) one gets 


T B(s)C({s) - A(s)D(s) 
S(s,x) = 8S, = x | P(s,x) Q(s8,x) | 


xXB(s 8,0 xA(s)D(s 


P(s,x)Q(s,x)  P€s,x)Q(s,x) 


ro P(s,x) = P(s,0 8,0) = 8X 
P(s,x)Q(8,x) 
Z es — om (4-52) 
Q(s,x) P(s,x) 
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Now if the variable x changes by an incremental amount x, then an 
approximation given by the first terms of a Taylor series, T(s,x), 


changes by an incremental amount given by 


AT(s,x) = T(s,X + AX) - T(s,x) 


4 AEG) bx 
& (42-) 1(2,x)5(s,x) (4-53) 


For a given input function R(s), the system normal response is 
given by ¥(s,x) = Rs) Tds,x) 
The change in responce, Y¥(s,x), caused by the increment Ax is 
given by AY(s,x) = Rs) AT(s,x) (4+54) 


Substituting Equation (4-53) for AT(s,x), Equation (4-54) becomes 





A¥(s,x) = (-4=-) R(s)T(s,x)S(s,x) (4-55) 


Equation (4-55) indicates that AY(s,x), the increment in the response, 
may be obtained by calculating the response due to an equivalent input 
perturbation signal, AR(s), applied to the original transfer function, 


T(s,x), where 





AR(s) = (22. R(s) S(8,x) (4,56) 


Simulation: Equation (4-55) may be represented in block diagram form 
as drawn in Figure 4-3, where S(s,x) has been represented according 
to (4-52). From Equation (4-44), polynomials Q(s,x) and Q(s,o0) can 


be written in expanded form as follows: 


Q(s,x) = (¢, + xd.) + (c) + xd) )8 + ereenenen- 
+ (c +xd_ 4) OP ‘ae (4-57) 
Q(s,0) = rs + Cc S + One neomasn of cen af ae (4-58) 
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Simulation of the transfer function Q(s,0)/Q(s,x) may be accomplished 
on an analogue computer as illustrated in Figure 4-4, A similar 
simulation is possible for the transfer function P(s,o)/P(s,x) so that 
interconnection according to Figure 4-4 gives an analogue computer 
method for obtaining the change in the response, AY(t,x), directly as 
a time function trace for any input signal. 


Sinusoidal Steady State. 


In the sinusoidal steady state, 3s = jW, 


T(jW,x) = T(W,x) eF°M4>®) (4-59) 
where, 
ol} ImT (jw, x) 
tan @ =» ReT(jQ,x) (4-60) 
and 
dnt (j0,x) = frft(@,x)| + JO(0,x) (4-61) 


From Equation (4-61) and the definition of S(jW¥,x) as given by 


Equation (4-49) it follows that: 


Re[S(jo,x)} =» —SytCax) (4-62) 


Im[S(jo,x)] = A8(2,22. (4-63) 


Using Equations (4-62) and (4-63), if then x changes by an 
increment, Ax, the corresponding incremental changes, Atu|T (4, x) | ; 


and A@(4),x), are given by 


a dx[t@,x)l - * 
(4% ) Re(S(j0,x)] (4-64) 
x 
40(O, x) 


taefay EMSC) ] (4-65) 
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Iv-4 - Sensitivity Integrals 

The expression for sensitivity leads to functions of parameter 
x (or k) and s, t or >. From an engineering point of view it would 
be convenient to have an overall measure of sensitivity which is 
independent of either time or frequency. Such a figure would be use- 
ful for comparing the relative sensitivities of different parameters 
in a system as well as presenting the possibility of including a 
sensitivity criteria for purposes of design or optimization. 
(a) - Phase Integral. 

Consider first the area of the equivalent input perturbation 
signal, Ar(t,x), relative to the area of the input signal, r(t). 
The total area of the input signal is given by the integral, 

: | 1 
I, = lim { r(t)dt Lim pa [Bir | 
t»O o C0 g 


= lim _sR(s) | R(0) (4-66) 


8S—>O 8 


In the above calculation use is made of the fact that, 


o 
and the final value theorem is applied. Similarly the total value of 


the area of the equivalent input perturbation signal is given by 


t 
Al, 2 lim 5 A r(t)de 
t 7c oOo 


ze As _R¢s)S(s8 5x) | 
ia SIGE) 
Ax 8 
= (4x) Lim er 6s)S (85x) 
$3090 


(42) R(o)S(0,x) (4-68) 
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If I, # 0, then the ratio of Equation (4-67) to Equation (4-66) 


denoted by I(x) is, 





AL 
I(x) s ey ' Ge 


) S(o,x) (4-68) 
i 


When I, = O it mesns R(0) = O and therefore from Equation (4-67) 
it follows that, Al, = 0, provided S(o,x) ¢@. When R(O) = 0 and 
S(o,x) =%, then At, may be evaluated by calculating the limit of 
the product R({s) S(s,x) as indicated by Equation (4-67). 

Alternately I(x) may be obtained by considering the area of the 


change in the response due to the increment 4x. Prior to the para- 


meter change the area of the response is given by 


I = lim y y(t,x)dt 


tom Oo 
= lim [% R(s) T(s,x 
t > co 8 
= lin s 1s) ¥%4s.x) 
s™o Ss 
= R(o) T(o,x) (4-69) 


After the increment x, the area of the change in the response is 


with the aid of Equation (4-54) given by 


t 
Lim § Ay(t,xddt 
tO Oo 


Al. 


(43) tin Ss RGEC 8 ex)S( 8.) 
x s 


(42) lim = R(s) T(s,x) S¢s,x) 


(42) R(0) T(o,x) S(o,x) (4-70) 


If I #0, the ratio of Equation (4-70) to Equation (4-69) exists and 
) 


is given by 
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AL | 
I(x) = = a (a) $(0,x) (4-71) 


© 





When I, = 0 it follows from Equation (4-80) that I, = 0, provided 
S(o,x) #0. When I, = 0 and S(o,x) scothen AL, may be evaluated by 
calculating the limit of the product R(s) T(s,x) S(s,x) as shown by 
Equation (4-70). 

Now consider simusoidal steady state amd the area of the 
incremental phase characteristic E(x.) plotted on a logarithmic 


frequency scale, defined by the integral 


@ 
E,@,x) = JA OW, x)d(e&) 


= d 
» fi e(,x) 2 (4-72) 
=a 


Substituting Equation (4-65) into Equation (4-72) we get 


aa da) 
1,@,%) = J) sx) Im(S(i0,*))] —3- = ’ 
eae : Ww 
x et SC52)} an (4-73) 


From the reactance integral theorem (Appendix IV), it is known that 


if H(ja) = RW) + jx) 
and lim H(jW) = O 
L) - CD 
then a2 
JS x(w) 2 | -4 lim H(s) (4-74) 
2% s=—O 


Applying Equation (4-74) to Equation (4-73) we get 


-4 (43) lim $(8,x) 


s 20 


wT [Ax 
- ¥ (42) 0,2 


= = (x) (4-75) 


a 


Ty (x) 


a 
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Equations (4-68) , (4-71), and (4-75) enable the following theorem 
to be stated. 
Theorem - For an incremental parameter change, Ax, the product of 
(4x) times the logarithmic sensitivity function evaluated at S = 0 
is oat to 
(a) - The ratio of the area of the equivalent input perturbation 
signal to the area of the input signal, provided the 
latter area is non-zero. 
(b) - The fractional change in the system’s response to any 
input signal, provided area of the response is non-zero. 
(c) = = * times the change in area of the sinusoidal steady 
state phase characteristic plotted on a logarithmic 
frequency scale. 
When S(o,x) = 0, Equation (4-75) indicates that the increment Ax, 
causes zero net change in the area of the sinusoidal steady state 
phase characteristic plotted versus dnt). Similarly, when S(o,x) = 0, 
Equation (4-70) indicates that 4x provides zero net change in the 
area of any time response characteristic which initially has finite 
none-zero aréa. 
If the area of response is zero and S(0,x) is finite there will 
also be zero net change im area produced byfx. S{o,x) may be cal- 


culated directly from Equation (4-52). 


~ 269,09) _ _F(o,9) 4-76a 
S(o,x) 0(0,x) P(o,x) \ ) 
<6 ee 
c + xd ” a +& xb (4-76b) 
re) oO oO ° 
where, Q(o,x) = cit xd | 
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P(o,x) = a: xb (4-76b) 


From Equation (4-76b) it is obvious that S(o,x) will be identically 


zero when a d zs cb, 
° Oo ° Oo 


(b) - The Magnitude Integral. 


It is also of interest to investigate the area of the change in 
the sinosoidal steady state magnitude characteristic as given by 


Equation (4-64). Consider the integral, 


oS 
AL, @,x) -) A twT(W,x) dw 


> 
= (42) § Re [S(jo,x)]dw (4-77) 
Oo 


From the resistance integral theorem (Appendix IV) it is shown that if 


H(j®) = R(j0) + jx(t) (4-78) 
and lim H(jo) = O 
W~» CO 
Then aad 
) R(©) dio = 5 lim sH(s) (4-79) 
° 80D 


Substituting Equation (4-79) into Equation (4-77) we get 





AL,(,x) = (42) 3 lim 8S(s,x) 


¥ 8 > co 
The use of this integral is illustrated in the following discussion. 
(c) - Sensitivity Squared Integral. 
In many cases S(o,x) is zero so it is not always a convenient 
criterion for comparing the sensitivities of different parameters in 
a system. This suggests the possibility of investigating the areas 


of the incremental error signals squared. 
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In previous sections of this chapter it was shown how the critere- 
ion of minimum square error to the incremental response was utilized. 
A sirpler and less restrictive approach is to investigate the area 
of the equivalent input perturbation signal squared, Thus, the 


following integral is defined, 


cS 


17 (x) = (or(t,x)] “de = [ar(t,x)] dt (4-81) 


re) o © 
where it is assumed that Ar(t,x) = 0 for t< 0. Applying Parceval's 


theorem to Equation (4-81), with the aid of Equation (4-56) we get 


I(x) = = [R(s)S(s,x)] [R(-s)S(-s,x)]ds (4-82) 
5 :) Baen 


when R(6) = 1, r(t) = O(t), then Equation (4-82) becomes 


2 fie 
17(x) = sx.) S(s,x)S(-s8,x)ds (4-83) 
3 5. 


Integrals of this type are tabluated in the literature (reference 2) 
and form the basis of analytical design procedures. In particular this 
integral can form a basis of comparison for the effects of different 
parameter variations upon system performance. The approach is ili- 
ustrated by a simple example. | 
Example: 
Consider a second order positional servomechanism subjected to 

a load disturbance torque T(t) applied to the torque shaft. The 
transfer function relating the output angie to the applied torque is 
given by ; 

or ~ Sees 7 
L Js + FS + K 


Where J, F and K are moment of inertia, viscous camping and gain 


10° 


Y¥(s,x) 


oY (8.x) 
P(s,0) 
P(s,x) | Ax/x 
~ AY(s,x) 
= 
at 
Q(s, 9) 
Q(s,x) 
¥Q(s.x) 


(\x/x) 


Fig. 4-3: Block Diagram for Generating the Incremental Response 
Caused by the Parameter Change, 4x. 
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Y¥(8,x) 


Fig. 4-4; 


Simulation of the Transfer Function 


Re. . (—E_) 
Q(s,x) \ Q(8,x) E 





A 


C1) | p— 


Y.(s,x) 
(4x/x) 
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respectively. It is desired to compare the sensitivity of the 
impulsive response when T, (t) 2 3(t) for incremental variations 


of F and K. From the polynomials of Equation (4-84), it follows: 


x s F 

P(s,F) = 1 

P(s,o) = 1 

Q(s,F) = x + Fs + K 


Q(s,o) = Cen K 
Substituting into Equation (4-52), we get 


Fs 


S(s,F) #= a (4-85) 
Js + Fs +K 
Similarly it can be shown that 
S(s,K) = -K (4-86) 
Je” + Fa + K 


Substituting Equation (4-85) and Equation (4-83) we get 





2 ar \? JF F 
nr) - (5 -Fs 8 de 
-j\° >Js + Fs +K Js ~- Fe+kK 
2 2 
a F AF AF by 
2J F ) a F ae 
and similarly 
2 2 
2 K /aK Wn ax.) , 
where 0 = K/3 as 
n 3 - F 
2V KJ 


Ww, is the undamped resonant frequency and S is the damping ratio. 
Equation (4-87) and (4-88) show that both sensitivity squared inte- 
grals are proportional to O- For 5> . a fractional variation of 


F is seen to have a greater effect upon its sensitivity squared 
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integral than the identical fractional variation of K ii:s upon its 


integral. 


IV-5 - Signal Flow Graph and Sensitivity 

The general formulation of the transmittance of a signal flow 
graph described by Mason, can be extended to evaluate the sensitivity 
function. Mason has derived a simple and very useful rule for eval- 
uating system transfer functions (or sometimes called transmittances) 
from signal flow graph. The transmission from any independent node 


x, called source,to a dependent node y, called sink, is 


SP AO 
" s xy xy 
xy A 


Where the system determinent A is 


Aw 1 -2L.¢ 2 Re - 2 bP viet + <--ereeg0s 
(ie, Oe i ahiaten * 


In Equation (4-90) L, represents a loop transmission, i.e. the 


transmission around any closed-loop and t Ly is the sum of al? loop 


transmissions. Li L! is the product of the loop transmissions of 


J 


any two loops that do not have a node or branch in common, and 


> L! L' 4s therefore the sum of all such products. Li L" Ly is 
im * J j 


the product of the loop transmissions of any three loops that do not 
have any branches or nodes in common. 


in Equation (4-89) is as follows: 


The definiton of P 4A 
xy xy 


Py is any direct transmission from the input x to the output y. Each 


Pins h 
Asy is associated with a particular Pay The value of oe is the same 
as A except fcr the removal of all terms containing any branches or 


~* "s,s set ay» a ii 
medes thet zpperr in FP. 
Nd 
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Fig. 4-5: 





A Signal Flow Graph for Applicaton of Mason's Gain Relation. 
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Example: 


Find A, Pio? as and Bio in Figure 4-5. The sum of the indivi- 


dual loop transmissions is 


Lb, = 1, 173 Tio + T3 7, Ts Ty, + 73 Ty Ts Te Tyo 


+ T, T, T, T, T, + T, T, T, (4-91) 


The above is best done in a systematic manner by starting with Ty 


and listing all loops containing Ty» then going on to To» and taking 


all loops containing T, that have not already been listed, etc. 


2 


Next we find ‘ae Li L. which is the sum of the products of the loop 
b 


transmissions taken two at a time, omitting all those that have any 


node or branch in common. To evaluate a L; ne the list in 3 L, 
5) & 


2 
is examined. The first term in i Ly is Tot 3t 19° The other terms 


are scanned to detect those that do not contain any of Tos Ty; or Ti0° 


The first term satisfying this condition is TyTetetzT > but Figure 


4-5 is examined to check whether these loops (ToT5T 16 and TT TTT 4) 


have any node in common. They do, and therefore their product is 


not included in L‘ L' However the last term T,T,T, has no branches 


i -j° 
in common with ToT 3T 16 and therefore Tota otstgrta appears in 
¢ ¢ 

Duy My Next we note second term in Lys Loe. , T7167), and examine 
as above. In this way it is found that, 

' ’ © 

2 Li Lt = (TjT Tyo) (TTIg) (4-92) 

Since pe Ly contains only one term, it is impossible for Z Lt Ly 


to exist. Then 


> _?f 2  o 
A 1 ci im Li by 


= 1 = (TTT yo + Ty%yTel yy + Ta2yTsTetig + TeTskeT7%o 


+ ToT eto] + ToT sl iotgtetg (4-93) 
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Continuing with evaluation of Pn and o we have, 


Fiom” Tit2T3*s*shen7 (4-94) 


Since each loop in tL, has at least one branch present in Pio? then 


A i= 1 (4-95) 


and Fe pa — T,T{T3T,T.T Ty (4-96) 


Substituting Equation (4-93) and Equation (4-89) we get Tyo" 


Formulation of Sensitivity Function: 





Sensitivity of a source to sink transmittance T,, with respect 


ij 
to a parameter is defined as 
5 13 = aa = oral (4-97) 
x Ux q x/x 
As discussed in chapter 3, the sensitivity function may be 


expressed in the following form: 








T T, ,(x=0) 
x F T F F 
x ij x x 


Where T, G9 is the transmittance from source i to sink j when the 
specified element x vanishes, i.e., when the branch x of the signal 
flow graph is removed and where Fe and a are the return difference 
and the null return difference respectively, both referred to the 
element x. Furthermore, 
2 yo: 
Th = — (4-99) 


as discussed in previous sections. Expressing Equation (4-99) in 
logarithmic form and differentiating with respect to bx, we get 


T Ww3(P. A,.) 
as iy bw : 


The terms on the right hand side are themselves sensitivity functions 
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and one may therefore regard the sensitivity function of Equation 


Yiw Sb 
lw x 


is the sensitivity of the graph determinent with respect to the 





(4-97) as a difference of simpler sensitivity functions. 


element, x. The way the determinent A was defined, it is then only 
a function of the feedback branch, taking the form of the sum of the 
products of branch transmittances,none of which appear more than 
once in any single product. Therefore we can write 4 as a linear 


function of s, 


and d 4 oe te 

a x 2S = ax = 4-J° (4-101) 
Where A° is the determinant when x = 0, i.e., when the specified 
branch x of the signal flow graph is removed, However accordinging 


to signal flow graph theory, 


F.. a nw (4-102) 


Dividing Equation (4-101) by A we get 





:- - 3 34/b_ = = ? = o l o 
ire 1- D/A as A 103) 


Similarly, the first term of the right hand side of Equation (4-100), 


which may be thought of as the path sensitivity, is found to be 








o 
cP, 4,.) 
<P | ie Ge her ta S ie ae (4-104) 
x D Py 5 fay? Fs 
Then Le 
S 4) sf = 54 = 3 zh e (4-105) 
x x x FE FS. 


Equation (4-105) was previously obtained in chapter 2. Also sensiti- 


vity could be written as, 
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A o 
5 13 a 4 2 (4-106) 
z Bay Mag 
In certain configurations this may be an easy way to evaluate 
sensitivity functions. However this is recommended as an interest- 


ing topic for further investigation to utilize the concept of signal 


flow graph and ite related techniques. 


IV-6 - Multiparameter Sensitivity 
So far we discussed the sensitivity of a transmission function, 
T, which depends upon a single parameter x and defined as 


st s d lyT 
x dhux 


(4-107) 
Now consider T as a function of n parameters, i.e., 


T ood T (x, » Xo» eevee egy x) 


Several definitons have been proposed. Hakimi and Cruz (see reference 
18) define multiparameter sensitivity by , 
maximum AT 
O< Ox, <% aT, J4 cilia 
where Q, is the maximum variation of the ith parameter as specified 
in particular application. For small J the fractional change in 


transmission is given by (neglecting second order effects), 


ae 5 om Te j 
7 d(T) L 5 (Lx,) d( lux, ) (4-109) 


Setting ws = nx. , we can consider the set of fractional parameter 


increments as a vector, 


and Equation (4-109) becomes the scalar product 
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Ste VY tht, --------- y}¢ dy (4-110) 


The sensitivity of transmission clearly depends upon the gradient 
vector VixT. Kuo and Goldstein (reference 18) suggest that multi- 


parameter sensitivity s" be defined as 


sy = J [mT(y,, ------ im) (4-111) 


This has the advantage that it agrees with the single parameter defin- 
ition in Equation (4-107). The magnitude of the multiparameter sensi- 
tivity function is 

| s? | 2 (Ste styl? (4-112) 


Donald A. Calohan (reference i7) defines the parameter sensitivity 





vector, 
A 
gt(3o) 50) dbwr ie (4-113) 
x \ kw X i 
i i 
7) eo ts pg VO ,s 
ma a i x, 
where 
A(wW) = bn. |r ¢5n) | , the attenuation 
B (d) = arg T(jW). the phase 
x, * element varying 
cS am a set of orthogonal unit vectors 
Hence, 
a A 
2 So xx, © Re (st) (4-114) 
i 4 
A 
5 ae Im (gh 65>) (4-115) 
x ii x 
£ i 
The magnitude of the vector re(st 3%), L.e., 
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7) 
Re Ga!» - 2 cos ” (4-116) 


is therefore an indication of the sensitivity. Multiparameter sensi- 
tivity could best be discussed in terms of state space, and is recom- 


mended as a topic for further work on this line. 
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V 
CONCLUSION 

The object of this thesis as stated in the "Introduction," was 
a literary research on the subject of sensitivity applied to the 
linear feedback control system, with only one parameter changing. 

Two different cases were discussed. One when the parameter changes 
are incremental and the other when parameter changes are moderate 
or large. In each case a system sensitivity was defined associated 
with methods of measuring this sensitivity function. 

In the case of incremental parameter variation, the concept of 
root sensitivity was presented and from that a graphical design tech- 
nique discussed. When the parameter variations were large, a new 
definition of sensitivity was introduced and it was shown that when 
the leakage transmission between the input and output is negligible 
the two definitions are identical. A design philosophy, and the use 
of Nichols Chart was also presented. Finally, the new aspects of 
sensitivity such as specification for sensitivity function based on 
forced response criterion and least mean square error were explained, 
A new approach based on the signal perturbation theory together with 


sensitivity integrals were also introduced, 
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Appendix I 

Bilinear Theorem, Return Difference, Null Return Difference 
1 - The Fundamental Feedback Flow Graph. 

The signal flow graph for a linear feedback control system is 
shown in Figure Al-l, where, 

tut is transmission from input I to output 0 

Le is transmission from source S to output 0 
ta is transmission from input I to control C 


tis is transmission from source S to control C. 


Then we can write, 


Output 0 = twit + tos? 
and S = ke ,/Q = kt.) 
Hence, 
kt .¢ 
0 D ci os 
7 T(s) tua a THE kt (Al-1) 


Equation (Al-1) is called the "fundamental feedback equation" and is 
the most important equation in linear feedback theory and all the 
significant properties of linear feedback systems are derivable from 
this equation. As an analytic tool it shows how a linear circuit prob- 
lem can be broken up into two simple circuits. 
2 - Return Difference. 

The system transfer function as derived above is 


Tey oe t _ + kt 


I oi ci°os ) 


/a - kt. 


suppose it is asked "what is the effect of the feedback around k, on 
the system function TT?" Clearly the effect is given by the quantity 
1 - kt, This quantity 1 - kt. has a very simple physical interpre- 


tation readily seen by opening the closed loop anywhere, for example 
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cs 
Fig. Al-l; Fundamental Feedback 
Flow Graph. 
t 
of 


cs 


Fig. Al-2: Physical interpretation of loop transmission function. 


t 
cs 





Fig. Al-3: Derivation of null return difference. 
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at AA' and measure the signal returned at A, the latter is kt. 
The difference between the injected and the return signal is pre- 
cisely 1 - kt. and is called the return difference for element k. 


yr 8 1 -kt 
Cc 


k (Al-2) 


8 
3 - Bilinear Theorem: 

The transfer function of any two port linear systems can be 
written in the form of the fundamental feedback equation (Al-1), such 


that all t, are independent of the system parameter k(s), providing 


j 
k can be represented as a controlled source S = kC and providing k 
appears only in S and nowhere else. 

The above theorem is called bilinear theorem because the funda- 
mental feedback equation may be considered as a bilinear transforma- 
tion of k. 

4&4 - Null Return Difference, Figure Al-2. 

The null return difference is the return difference evaluated 
under the condition that the input is adjusted to give zero output. 

In terms of the signal flow diagram the null return difference 
is determined as follows. The diagram is broken at A and A‘ and a 
unit signal is transmitted from A'. The signal reaching S is k. This 
is transmitted back to C and also to the output 0. If the input is 
adjusted to make the output zero the signal “kt. must be arriving 
at O along the branch from I. Hence I is okt /tos: This input 
is also transmitted along the branch from I to C and the return signal 
is the sum of the two signals arriving at C. Then the null return 


difference with reference to k is given by the relation, 


' BS @ eS 
Lat : (kt. kt toy /toy) 
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The system transfer function can be written in terms of return diff- 


erence and null return difference, 


kt .¢t Tl =kt) . # Re t-te 
cs 


ci os cioos oli 
T ©) t., + So 
oi 1 - kt oi l - kts (A1-3) 


Where Fy is the null return difference and Fe is the return diff- 


erence with reference to k. 
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Appendix II 
The Curvature of the Root Locus at Ordinary Points 
The curvature of a curve, which is defined by means of an 


implicit relationship F(x,y) = 0 is given by 

















2 2 
ji. *, ee A ae eee Lor) 
Sf Fr? A aes 
where F.. a oa (A2=2) 
r 
F, * (A2-3) 
2 
UF \“F 
Pex s ~- (22). 852 (A2-4) 
2 
eae Ck ae: SF \_ °F 
Fay dy Vx 3x sy y xdy (A2~-5) 
2 
F MF 
Py ay (2z-) ay? sie 


If we introduce a new cartisian system with its origin at s, and 


j 


oriented so that the tangent to the locus coincides with the u-axis. 


Pe 


Then in this system for the implicit function F(u,v) we have 


(Figure A2-1) 


SF 
TG = F (0,0) = 0 
v20 
u=( 
and 9 J 
l F F - 2F FF +F F 
ae uu Vv uv uv VV U (A2-7) 
Pp aye + yp? 


and since F (0,0) = 0 and curvature is evaluated at point (0,0), then 


1 . uu 


af — us(Q 
v=0 
The center of curvature is on the line normal to the curve. Because 


of the perpendicular choice of coordinate the positive sign of / means 
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S-Plane 





Fig. A2-1: UU, Vv coordinate system. 





Fig. A2-2: (S - oy in v, and u coordinate. 


125 


that the center of curvature lies on the positive v-axis and vice 
versa. Since the gain, K(s} is a rational function of s and analytic 


except at its poles s = Zs then 


n m 
Ky(s) = Av (K(s)] = Di bw(s-P,) - 2, dw(s-2,)+ 05 
j=l i=] (A2-9) 


This is a two dimensional complex potential function such as would 
be used to describe the field around a set of n unit positive line 
charges at the points Le and m unit negative line charges at the 
points Z, for such configurations, the lines 

Re K, (s) = An| Ks) | = constant (A2-10) 
would be equipotentials, and the lines 

Im K, (s) = argK(s) = constant (A2-11) 
would be stream-lines or lines of force. Since the above relations 


and the relation 


m n 
arg K(s) = v2 arg (s - Z,) - 2 arg (s - P,) 
j=] j=l (A2-12) 


are equivalent with a particular value for the constant, it follows 
that the root loci coincide with certain lines of force in potential 
analogy. 

Having the above idea in mind we can define function F as F = 
ImK and sometimes is defined as F = Im Wx, which is more suitable 
for our purpose. As K takes on positive real values, F is equal to 


zero and therefore serves as the real function of Equation (A2-1). 





Therefore ‘ F (Im 4K) 
- - 2g tp 
P oe (1m 4K) _ 
In [« nk) uu 
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In general K is expressed as a function of the variable Z = u+ jv 


and since 1nK is an analytic function then, 








+ * Jarier aia? 
Hence 
1 Im( JK) _ Im(4» K)_, 


% " "WmiGK, ~ ~ Ree, oer) 


To evaluate Equation (A2-15) we assume Z = (8 « s ye i%, and for any 


function £ we have 


























df df, ds © ae ay 
dZ ds dZ ds 
jae coe dae (i) , gf _d%s 
azz d 2 dZ ds dz- 
d*£ 25% 
= e (A2-16) 
2 
ds 
and similarly 
d'£ d"¢ (-)* vf jn? (A2-17) 
az” de" dZ ds 
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Appendix III 

Derivation of an Expression for the Sensitivity of a Multiple Order 
System Root 

If we write the expansion of the total differential of L(s) to 


include higher-order terms (see reference 15, page 81), we have 





(dL), = 0 ® $e do+ Sh ax +2 Ub az + FL ap 








S IK i 8Z, i 
os a ‘ 
x same as above bracket| 
iM 5 [ ai ] oe I 
where oe 
SL SL . >» *y a 
= ds + yK dK are (ds) 
ne? a-l a(a-1) eae a-2 2 
+ a-——— (ds) aqK + - ———_ (ds) (dK) 
¥ 82 yx Sa ili’ 
fe oe (ax)? 
axe 


Next retain only the lowest order terms for each parameter and 
note that the first (N-1) derivatives of L with respect to S are zero 


at S = S.,. Then 





j > 
) 1/N 
d5 etter ig: SL ax + Z 2h az, +7 2 aP, 

J _ 3 1 2% i °"4 
8=8 
j 

This suggests the following notation: 
1/N 


8s 8 
ds. = Wee; “* § 2 asc® Gees’ ap 


and comparing the above two equation we have 
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Appendix IV 
Resistance and Reactance Integrals 
To prove the resistance and reactance integrals, we start with 
Cauchy's residue theorem. Consider any region defined by a closed 
boundary in which and on whose boundary H(s) is single valued and 
analytic, except for poles of any finite numbers of multiplicity. 
Then the line integral of H(s) around this boundary is related to 
the residues of those poles of H(s) which are located inside the 
region as follows, 
& H(s)ds = 2113 2 residues. 
If H(s) has no singularities in the right half plane and as on 


Figure A4-1, then 


jR -V/2 
& H(s) ds = lim J H(s)d(jw) + lim f H(R,0)d(Re?”) 
R20 -jR Reo 11/2 
= I, + 1, = 0 (A4-1) 


Next we consider some properties of the loop transmission function 
H(s). The loop transmission function can be regarded as transforms 


of impulse responses of real systems. The transform of a time function 


is a kt 
H(s) = h(t)e dt (A4-2) 
So 
oO 
H(ja) = if h(t) [cosHt - j sinYt]) dt 
= 
= R(wW) + jx (a) (A4-3) 
where sal 
R(W) = h(t) cosWt dt is an even function of W 
a 
and x@) = f h(t) sind&t dt is an odd function of ®. 
= 00 
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S-Plane jR 


-jR 


Fig. A4-1: Application of Cauchy's residue theorem. 





Fig. A4-2: Contour for F(s)/s. 
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Thus, H(jo) = RW) + Jx(W) = H(-jw) = A(-) - Jx(-w) (A4~4) 
It is often necessary to focus attention on H(s) for s near zero and 
for s near infinity. Expressing H(s) about the origin gives 


= t t 2 | 3 aanauawea ‘— = 
H(s) A, + B's + A 53 “+b B" 48 + (A4-5) 


where all the A‘, 8 and B’ 8 are real. The notation for the expansion 
of H(s) at infinity is 


A). seal ee (BESS) hee (AM Jan) deuoaase-== (A4-6) 


where all Ay, and Boo are real. 

Consider H(s) so that it satisfies Equation (A4-4), is finite 
(or zero) at infinity and is analytic in the right half plane. Apply- 
ing Cauchy's theorem to [H(s) - H ()] over the right half plane 


whose boundary is the jw-axis and the right half-infinite semi-circle, 


) [H(s) - H(W)jJds = O 


R 
= lim ) [H(J0) - A] jd 
R2© -R 
(1/2 
+ lim (H(Re!®) ~ a,j d(Re!”) 
Ra -1/2 
(A4-7) 
But let H(jwW) = RW) + jx(W&) , and note that, 


R(®) = R(-6) 
x() = -x(-W) 


Consequently, 


R (sa) 
lim s [HC j0) 2 A.J jdw «a 12) [R@) - AL) dt 
R?@ -R 0 


As s approaches infinity the only term of [H(s) - A,] that contributes 


to the second integral in (A4-7) is (By,/s) and then 


» ia -1i/2 
lim [H(s) - AJds = lim § BL /e ds 
RPO -71/2 R> 7/2 
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= 11/2 


= lim JS Bin /s Re!® de 
Re 71/2 
fi 
= lim Bi, jd9 = #£~ iB? 
Rao TI/2 a - 


and therefore, 


j2 Sra - a2 dw - jmBL, = 0 
0 


or a 


J [R(O) - AJ dw = 1/2nB, (A4-8) 


° 
Note, that if we had used H(s) rather than [ncs) 2 Add, as the integrand, 
then the integral will be infinite unless A,j= 0. If 


lin H(s) = A_ = O 


CO 
Sa 
then Bo = lim SH(s) 
8 > co 
Ce 
and 
R(W) dW = 1/2T1lim SH(s) (A4-9) 
0 s> co 


This is the simple form of “Resistance Integral" used in Chapter IV. 


Reactance Integral: 


We apply Cauchy's theorem to H(s)/S over the contour Figure A4-2 





R et/2 A 
& 1g 5 aH 4. Wf iG a(re?®) 
: r»oO /f ) n/2 Re 
R-» co 


r “2 A 
+ J FGM) aw+ f Oe ggeel"S Caggesion 
@-R LJ -1/2 ne 


When the first and the third term of the above equation are combined, 
the odd parts cancel and the even parts add. Writing H(ja) = R(% 


+ jx(W), Equation (A4-10) reduces to 
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2 x(o) $2 + m(A- A) = 0 
#- ©O 
or > 
¥ x(W) “2 = -T7/2 (A, - A) (A4-11) 
= 0 


If A =# QO, then Equation (A4-11) reduces to 


fe a) 
x) a = -W/2 lim H(s) 
= 0 s30 


which is the simple form of reactance integral used in Chapter IV. 
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